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Preface 

This thesis is devoted to the investigations of gravitational wave data analysis 
from a continuous source e.g. a pulsar, a binary star system. The research 
in the field of gravitational radiation started after its formulation by Einstein 
in (1916) as propagating gravitational disturbance described by the linearised 
field limit of his General Theory of Relativity (GR) but has received seri- 
ous impetus towards its detection after the announcement of its detection by 
Weber in 1969 using aluminium bar detectors. This field has now emerged 
and established itself with General Relativity (GR), Astrophysics and Nu- 
merical analysis as its equally important facets. Ofcourse, the technological 
advancements being employed in the construction of detectors with day by 
day improving sensitivity have played the crucial role. To date, although the 
results of Weber could not be confirmed and we do not have as yet any direct 
detection of gravitational wave (GW) yet it is not a matter of concern. Be- 
cause on one hand the sensitivity required for the announcement of definite 
detection of GW bathing earth is yet to be achieved by the detectors whereas 
on the other hand we have an indirect evidence of the existence of GW ob- 
served in 1974 as the slowing down of the binary pulsar PSR 1913 + 16 arising 
because of back reaction of GW emission. GW Scientists all over the globe are 
putting more persuasive arguments regarding the feasibility of GW detection 
in 'near future' and the advantages to be achieved once the "Gravitational 
Wave Astronomy" as they call it, is established. A huge amount of money is 
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involved in these projects to the extent that many of the detectors are built 
in collaboration e.g. American Laser Interferometric Gravitational Wave Ob- 
servatory (LIGO), Italian French Gravitational Wave Observatory (VIRGO), 
British-German Observatory (GE0600). As a consequence the literature is full 
of update reviews on gravitational wave astronomy notably by Thorne (1987), 
Blair (1991), Schutz (1999), Grishchuk et.al (2000), where the related issues 
viz., the fall outs, pre- requisites and prospects are discussed and scrutinised 
with rigour and minute detail. 

The thesis starts with two introductory Chapters dealing with gravitational 
wave and data analysis concepts, respectively. The matter covered in these 
Chapters are restricted to the extent they are supposed to provide continu- 
ity and coherence to the investigations presented in Chapters 3-5. Further, 
the material presented in these chapters depend heavily on the reviews by 
Grishchuk (2000), Schutz (1999), Sathyaprakash (1999). Our job have been, 
frankly speaking, to cut and paste. Let us hope this has not defaced the truth. 
For the portion describing GR we have adopted notations and conventions of 
Schutz (1989). The source codes of the numerical computations are relegated 
to the appendices A-C. 

The term "Pulsar" used in the thesis refers to spinning neutron stars. A 
Pulsar will emit GW signal over extended period of time only when it has long- 
lived asymmetry. Several mechanisms have been given for such an asymmetry 
to arise. Some Pulsars emit GW almost monochromatically and are remark- 
ably good clocks as its periods have been measured upto 13 significant digits. 
However, the GW signals from Pulsars are very weak (< 10~ 25 ) and will be 
buried in the broadband noise of the detector. In order to detect the signal 
from the dominant noise one has to analyse the long time observation data 
of months/years. The output of detector is highly involved function of many 
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initial parameters. It is usually not possible to obtain the Fourier Transform 
(FT) analytically. Hence, FT has to be obtained via numerical methods. But 
it appears to be computationally demanding even for the standard computers 
expected to be available in near future. Therefore, one will have to work with 
efficient data analysis techniques, efficient not only in picking weak signal from 
the noisy data but also in terms of computing-cost requirements. Chapter 2 is 
on the problem and technique for the data analysis. The noises and sensitivity 
of the detector has been briefly described. The Matched Filtering, a technique 
of the optimal method for detecting unknown signal and which describes drop 
in signal-to-noise ratio in terms of Fitting Factor (FF), is discussed. The 
computational cost and detection criteria are also discussed. 

The output of detector will be Doppler modulated (both in frequency and in 
amplitude) due to the motions of Earth and star. To work in Fourier space, one 
usually employs the numerical method called Fast Fourier Transform (FFT). 
Fortunately, we have been able to develop analytical FT. We have employed the 
results to obtain efficient data analysis algorithms for continuous gravitational 
wave (CGW). 

In Chapter 3 we have studied the noise free complete response of Laser 
Interferometer detector for CGW for its arbitrary location and that of the 
source. We have taken into account the rotational motion of Earth about 
its spin axis as well as its orbital motion around Sun. We have developed 
analytical FT of Frequency modulated output for one day observation time. 
In Chapter 4 this analysis has been generalised for (i) one year observation time 
and for (ii) any arbitrary duration of observation data set. The emission of GW 
from pulsars as such may contain two or more frequencies. Hence, we finally 
generalise the transform for N-component signal. We also give the method to 
account for spin down of CGW. The Frequency modulation (FM) smears out 
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a monochromatic signal into a small bandwidth around the signal frequency. 
The FT contains double series of Bessel functions and the computational cost 
depends primarily because of Bessel functions which in turn depends on source 
location and frequency. 

The Amplitude modulation (AM) arises due to the anisotropic response of 
the detector i.e. the detector possesses a quadrupole antenna pattern. The 
effect of AM is to split the signal frequency f a into five lines corresponding 
to frequencies f Q ± 2f rot , f ± f rot and f a where f rot represents the rotational 
frequency of Earth. Hence, for optimal detection of Amplitude modulated 
signal one will need four additional linear filters. For a 2-component pulsar 
signal the number of filters will increase further by a factor of 2. The cases of 
Doppler modulation investigated by us reveal that as a result of the AM the 
most of the power lies in the frequency band / + 2f rot and the least in the 
frequency / - 2f rot . 

An important advantage of our analytical results over FFT arises when 
one looks for the resolution of the FT. The resolution offered by the FFT is 
inversely proportional to the observation time, T Q whereas, the analytical FT 
is independent of T . This means that for short observation time the resolu- 
tion offered by FFT is small whereas for long observation time there would be 
high computational demands. We have developed, using our analytical FT, a 
semi-analytical formula to calculate the Doppler shift for one year observation 
data for any source location and frequency. A monochromatic signal of one 
kHz frequency from a source located at 6 = n/2 is found to spread into ap- 
proximately 3.13 millions Fourier bins. We observe that the number of Fourier 
bins increases with the source location 6 and frequency f Q . 

The strategy for the detection of GW signal is to make use of FT to dig out 
the signal from the noisy output of detectors. To achieve this, one constructs 
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the templates which are best educated guesses of the expected signal wave- 
form. In Chapter 5 we have applied our results on FT obtained in chapter 4. 
Applying the technique of Matched Filtering we have computed the number 
of templates required for all sky search problem for a pre- assigned (FF). It is 
estimated that one would require about 1.44 x 10 10 , 3.5 x 10 10 and 5.5 x 10 10 
number of templates for respective data analysis of 30, 120 and 365 days for 
a FF of 0.97. We have studied the variation of the estimated number of tem- 
plates with frequency. It is found that the numbers of templates required 
are approximately 1.22 x 10 10 , 2.16 x 10 10 and 5 x 10 10 for respective signal 
frequencies 20, 50 and 100 Hz for analysing data of 120 days with a FF of 
0.97. In this Chapter we have also studied the variation of the FF with source 
location [9, <p} and have found a marked symmetry. The FF for any arbitrarily 
chosen 9, is the same as corresponding to n — 9. Similar symmetry have also 
been observed for 0. This symmetry property will reduce the computation 
burden by a factor of four for all sky search problem. It is not clear whether 
this symmetry property can be established analytically as well. However the 
source location, because of this symmetry, becomes uncertain and its exact 
location determination will require t he use of other methods. 



The thesis has been typeset by the candidate in BTgX. 
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Chapter 1 

Introduction: Gravitational 
Waves, Detectors and Sources 

Gravitational waves (GW) like any other type of waves are propagating per- 
turbations of some flat background space-time. These are identified as small 
ripples rolling across space-time in the same manner as water waves are on an 
otherwise flat ocean. These waves originate from the most energetic events in 
our Universe such as rotating neutron stars, colliding neutron star binaries, 
supernovae explosions and gravitational collapse in Black Holes. They mani- 
fest themselves as strains in space-time that periodically stretch and compress 
matter. The GW emanating from a binary may be represented as in Fig. 

Newtonian gravitation does not have the provision for GW. In several 
Lorentz covariant gravitational theories e.g. scalar, vector, tensor theories 
the GW arise as the spreading out gravitational influence. The basis of most 
current thinking of GW is Einstein's theory of General Relativity (GR). In 
fact, GW was first studied by Einstein in 1916 by applying linear field ap- 
proximation to GR. However, he was misled to the result that an accelerated 
spherical mass would emit GW. He corrected his mistake in 1918 and showed 
that the first order term of GW was quadrupolar. We discuss the properties 
of GW in the next section. 
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Figure 1.1: Gravitational wave from binary. 



The strength of GW is so small that there is no hope of its detection in 
the manner Hertz demonstrated the existence of electromagnetic waves. The 
reason may be attributed to the extremely small value of the universal con- 
stant of gravitation (G = 6.7 x 10 _11 Nm 2 kg -2 ). The dimensionless amplitude 
of typical GW reaching the earth is only of the order of 10~ 17 . This put in 
other words means that such a GW will put a rod of one meter length into 
a oscillation with an amplitude of one millionth of the radius of a hydrogen 
atom. Any chance to observe the effects of GW would require acceleration 
of astrophysical size masses at relativistic speeds. Several GW-detectors are 
being built at present, or have been proposed for the future, that hope to 
detect these very small vibrations. The best known of these are the Laser In- 
terferometer Ground Observatory; LIGO, and the Laser Interferometer Space 
Antenna] LISA. As the names suggest, the former is a terrestrial observatory 
and the latter a space based one. Needless to say, what enormous engineer- 
ing achievements these detectors require to detect the oscillations mentioned 
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above. 

The interaction of GW with matter is also very small and this leads to 
some profound benefits for astrophysics. The wave will not scatter, so they 
emanate undisturbed from the regions of their origin viz., densest regions 
of the Universe, inner cores of imploding stars, the earliest instants of the 
primeval Universe, and the formation of Black Holes. Hence, they will provide 
information on the dynamics in these regions about which there are no other 
means. 

A piece of indirect evidence affirming the existence of GW came in 1974 
when Joseph Taylor and Russell Hulse studied the pulsar PSR 1913+16 which 
is a binary partner of compact neutron star and is orbiting around the other in 
a slowly decreasing eight hour period. The orbital frequency acceleration was 
observed to be in perfect agreement with what expected from the energy loss 
arising due to their calculated GW emission (Damour and Taylor, 1992). For 
their accomplishment, Hulse and Taylor were awarded the 1993 Nobel Prize. 

Unfortunately PSR 1913 + 16 is presently emitting GW too feeble for 
their direct detection on Earth and such a condition will continue until the 
end of its pas-de-deux inspiral some 350 million years hence. Interestingly, the 
density of observed binaries in our galaxy is such that a few late stage inspiral 
events are expected per century to occur. Hence, the study of the orbits of the 
binary neutron star system and their coalescence have played important role in 
the setting design criteria for some of the instruments to come into operation 
in the next few years. To the date we have only the indirect evidence of 
GW. Yet this has not been of much concern for theoretical physicists because 
the strength of the postulated GW signals are below the detection threshold 
currently available. 

The quest to detect GW started in its earnest as early as in 1965 with 
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the pioneering work of J. Weber on resonating bars. These are basically high 
quality bells designed to be rung by transient GW. With this beginning a 
small GW detection community has arisen and thrived continuing to improve 
on the original idea. The main effort is directed towards noise reduction by 
introducing ever lower cryogenic temperature and ever more sensitive displace- 
ment sensor to improve detector sensitivity. Many more bars; progressively 
sophisticated ones have been built. The Explorer which is quietly functioning 
at CERN since 1989, is one such detector worth mention. We discuss GW 
detectors in section 2. 

There is an excellent prospect of detection of GW in the near future. There 
will emerge what has been called "Gravitational Wave Astronomy" . This will 
provide another window for observing Universe. The expectation is that it will 
uncover new phenomena as well as add new insights into phenomena already 
observed in electromagnetic part of the spectrum. GW emanate due to the 
accelerated motions of mass in the interior of objects. As remarked earlier, 
these regions are otherwise obscure in electromagnetic, and possibly, even in 
neutrino astronomy. The GW arise due to the coherent effects of masses 
moving together rather than individual motions of smaller constituents such 
as atoms or charged particles which generate the electromagnetic radiation. 



1.1 Gravitational waves 



The gravitational waves like electromagnetic waves may be defined as time 
varying gravitational fields in absence of its sources. The gravitation theory, 
widely believed, is the one proposed by Einstein and famous as GR. In its 
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mathematical essence it is expressed via the equation known by his name viz. 

G Q /3 = Ra/3 — -^9af3^ = SvrToJll (1.1) 

The conceptual meaning of this equation is that a localized density distribu- 
tion characterised by the energy momentum tensor, T a p curves the space-time 
around it. The curvature and other geometrical properties of the space-time 
are characterised by Ricci curvature tensor R a p and the underlying metric 
tensor g a p. The curved space-time forces a free mass particle or light to follow 
the geodesies. These geodesies measure the effect of gravitation. The GR in 
its weak field limit yields the Newton's theory of gravitation. In this approxi- 
mation the gravitational field is considered to be represented by a metric with 
a linear modification of the background Lorentz space-time metric, r] al 3. It 
is important to remark that many of the basic concepts of GW theory can 
be understood in this approximation and were introduced and developed by 
Einstein himself. At present G W theory; its issues and problems, are dealt 
in detail in research text books of GR. Due to the limited scope of the the- 
sis we base our discussion on linear field approximation and incorporate few 
aspects of GW required for understanding of the subject matter presented in 
the thesis. 

1.1.1 Linear field approximation 

The metric of the space-time may be expressed as 

9a/3 = Va/3 + h aj 3 , \h aj3 \ < 1 (1.2) 

It is straight forward to compute Einstein tensor, G a p and one obtains 

1 For notations, conventions and basic definitions the reader may refer to Schutz (1989) 
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where 



h »P = h *P _ hq^h ; h = h a a = -h a a (1.4) 



We raise and lower the tensor indices using the Lorentz metric. Let us note 
that the expression for G a p would simplify considerably if we could require 

= (1.5) 

In fact, there is a gauge freedom available as infinitesimal coordinate transfor- 
mation defined as 

x a > X ' a = X a + £ a (x 13 ) , |f a |«l (1.6) 

where £ a is arbitrary. This gauge transformation preserves (|1.2|) . It can be 
shown that, in order to achieve the condition (jl.5|) . £ a has to be chosen as to 
satisfy 

where the symbol □ is used for the four dimensional Laplacian: 

d 2 



°f = f' t , ll = n^f^=\ < - w + v A )f (i-8) 

The gauge condition defined via (jl.5j) is called Lorentz gauge. We adopt 

this gauge. In literature harmonic gauge and de Donder gauge are other names 
for this gauge. The Einstein tensor, to the first order in h a p becomes 

G a ? = - l -Uh^ (1.9) 

Thus the weak-field Einstein equations take the form 

UhT = -16ttT^ (1.10) 
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In the Newtonian limit where the gravitational fields are too weak to produce 
velocities near the speed of light, 

|T°°| > |T oi | > \T j \ ; T°°~q, O 2 ~ v 2 (1-11) 

Equation (jl.lOj) now results into 

Xj 2 h°° = -lQng (1.12) 

This equation may be compared to Newtonian equation for gravitational po- 
tential, if i.e. 

vV = -4tt£ (1.13) 
One obtains, after some calculations 

h°° = 2<p , h xx = h yy = h zz = -2tp (1.14) 

and accordingly the space-time metric for a Newtonian gravitational field is 
represented via 

ds 2 = -(1 + 2tp)dt 2 + (1 - 2tp)(dx 2 + dy 2 + dz 2 ) (1.15) 

1.1.2 Propagation of gravitational waves 

The GW in the regions far away from its sources may be described in the 
weak field approximation. The Einstein equations, in absence of matter and 
electromagnetic fields {T a p = 0) take the form 

Uh ^= ( _ ^ + V a )^ = ° (1.16) 
This is the three dimensional wave equation and its solution may be taken as 

h a/3 = k^expiik^) (1.17) 
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where {A°^} are constant components of some tensor and {k a } are the con- 
stant components of some one form satisfying the equation: 

Tcf^k^kv = k v k u = (1.18) 

This means that 1)1.17)1 represents a solution of (J1.16)) provided fc 7 is null one 
form and the associated four vector k a is null. The value of h a ^ is constant 
on a hyper surface on which k a x a is constant i.e. 

k a x a = k a t + k.x = const. ; k={k i } (1.19) 

It is conventional to refer k° as w, which is called the frequency of wave: 

k a : (w,k) (1.20) 

The gauge condition (jl.5|) now requires 

A Q % = (1.21) 

which means that A a/3 must be orthogonal to kp. 

The solution ()1.17|) represents a plane wave propagating with velocity of 
light. In physical applications one has to consider the real part of the solution. 
Having obtained the solution h a/3 one still has the freedom of choosing specific 
£ a with the requirement that it represents some solution of the Eq. (jl.7j) : which 
in view of (jl.lfij) becomes 

("|r + V 2 )r = (1.22) 

Let us take its solution as 

e a = B Q exp(^) (1.23) 

where B a is constant. It can be shown that the freedom available in choosing 
the values of B a may be employed such that new A a p satisfy the conditions 

A a a = (1.24) 
A^U^ = (1.25) 



8 



where U a is some fixed four velocity. Eqs. (jl.21|) . (jl.24|) and (jl.25|) are called 
the transverse traceless (TT) gauge conditions. We choose U a as time basis 
vector U a = 5 a . Let us take the direction of propagation of the wave as 
the z axis of the coordinate system. Now using (jl.l8|) and (jl.2U|) we have 
: (w, 0, 0, w). Now Eq. (OS)) in view of (fHTTl) implies (i) A a0 = and (ii) 
A az = for all a. This is the reason the gauge is called 'transverse' gauge. 
Further the trace free condition (jl.24j) require A xx = —A yy . Now the non 
vanishing component of A a p may be expressed as 



( \ 

A+ A x 

A x -A+ 

\ ) 



A-xx — A + , A xy — A x (1.26) 



Thus there are only two independent components of A a p; A + and A x . Let us 



note that traceless condition (11.2411 results into 



h% = h T J (1-27) 

We have considered the plane wave solution of Eq We know that 

any solution of Eqs. p.21j) and f)l. 16|) may be expressed, because of theorems 
on Fourier analysis, as a superposition of plane waves. Hence if we consider 
the waves propagating along z axis we can put all such planes waves in the 
form (jl.27j) . Thus any wave has only two independent components ti^J and 
h X y represented respectively as h + and h x corresponding to A + and A x . 



1.1.3 The effect of waves on free particles and its po- 
larisation 

The independent components h + and h x represent the polarisation states of 
the wave. To understand their nature in little detail it is instructive to discuss 



9 



the effect of GW as it hits a free particle. Let us choose a background Lorentz 
frame in which the particle is initially at rest. We may employ the initial four 
velocity of the particle (U a = 6 a ) to define the TT gauge of the wave. A free 
particle obeys the geodesic equation 

d\] a 

+ r^U^U" = (1.28) 



dr 



This equation may be used to obtain the initial acceleration of the particle: 

fd\1 a \ 1 

= = (Vo + V,o - hoo,p) (1.29) 

In view of Eqs. ()1.17j) . (jl.26|) and (jl.27|) . the initial acceleration of the 
particle vanishes. This means that the particle will be at rest a moment later 
and consequently will be there forever. What does it mean ? Is there no effect 
of the GW on a free particle ? No, the interpretation presented at its face 
value is quite misleading. The result only means that the choice of the TT 
gauge employed resulted into a coordinate system for the wave which stays 
attached to the individual particles. 

To get a better measure of the effect of the wave let us consider two nearby 
particles situated at the origin (0, 0, 0) and on the x axis, (e, 0, 0) separated 
by a distance e. In view of our above discussion the particles remain at their 
initial coordinate positions. The proper distance between them is 

A I = J \ds 2 \ 1/2 = J \g aP dx a dx p \ l/2 
\g xx \ 1/2 dx « \g xx (x = 0)| 1/2 e 



Al = jl + i/^J(x = 0)}e (1.30) 

Thus the proper distance between two particles (as opposed to their coordinate 
distance) does change with time. The effects of the wave may also be described 
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in terms of the geodesic deviation of the separation vector, r] a , connecting the 
two particles. It obeys the equation 

^rf = RV^UV (1.31) 

It can be shown that for the particles initially having the separation vetor, 
r] a — > (0, e, 0, 0) one would get 

B 2 1 B 2 B 2 1 B 2 

—n x - -e—h TT —n y - -e—h TT (I 321 

3P V ~2 M* *" dP V ~ 2 Bt 2 xy 1 } 

Similarly, for initial separation vector, r] a — > (0, 0, e, 0) we would get 



B 2 1 B 2 1 B 2 

^_ n y = -e—h TT = —e—h TT 
Bt 2 ' 2 Bt 2 yy 2 • Bt 2 xx 



Let us note that in view of the results of the previous section we may write ()1.17j) 

as 

h a p = A a/3 exp(wt — kz) (1-34) 

Thus the separation vector rj a of the particles oscillates. 

Let us consider a ring of particles initially at rest in the x-y plane as shown 
in Fig. (|1.2l -a). Suppose a wave having h^J ^ 0, h^y = hits the system of 
particles. The particles will be moved (in terms of proper distance relative 
to the one in the centre) in the way shown in Fig. (|1.2l b). Similarly, for a 
wave with = = h^y, ti^y ^ the picture would distort as in Fig. (jl.21 
c). Since h^J and ti^y are independent, Fig. (|1.2l -b) and ()1.2< -c) provide the 
pictorial representation of the polarisation states of the wave. Let us note that 
the two polarisation states are simply rotated by 45° relative to each other. 
This is in contrast to electromagnetic waves where the two polarisation states 
are at 90° to each other. 



11 



(a) 



♦ 

(b) 



(c) 



Figure 1.2: (a) A circle of free particles before a wave travelling in the z 
direction reaches them, (b) and (c) Distortions of the circle by a wave with 
the '+' and ! x' polarisation. The two pictures represent the same wave at 
phases separated by 180°. Particles are positioned according to their proper 
distances from one another. 

1.1.4 Generation of gravitational wave 

To understand the generation of GW it is sufficient to discuss the weak field 
limit equation (jl.K)j) . rewritten as 



Let us assume, for sake of simplicity, the time dependence of T^ as sinusoidal 
oscillation of angular frequency u: 




(1.35) 



T 



(1.36) 



and look for a solution for h ai/ of the form 



h 



lip 



B fJ , v (x 3 )e 



(1.37) 



This, in view of Eqs. (jl.35|) and (jl.36J) . would require B^ v to satisfy 



( V 2 + ^ 2 ) B. 



16ttS 



(1.38) 
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Outside the source i.e where = we want a solution which would represent 
outgoing radiation far away. Hence we may take the solution as 



B M „ = [-^j e wr (1.39) 

Obviously, A^ u is to be related to S^. Under the assumption that the region 
of the space in which S^ v ^ is small compared to 2ir/u it can be deduced 
that 

This assumption is referred to as slow motion approximation since it implies 
that typical velocity inside the source, which is uj times the size of that region, 
should be much less than 1. All, except the most powerful sources, satisfy this 
assumption. Thus we get 

V = (£) (1-41) 

This means that the generated GW has frequency uj. This relation may be 
expressed in terms of useful quantities with the help of the following results. 

(i) The energy momentum satisfies the conservation equation 

T M % = (1.42) 

and obeys the identity: 

(ii) 



dt 2 



J T°°x l x m d 3 x = 2 J T lm d 3 x (1.43) 
(hi) The quadrupole moment tensor I lm is defined as 
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and which, in view of (jl.36J) . may be expressed as 

where D lm represents the time independent factor of I lm (Misner et al., 
1973). 

The first result gives 

J"° = => V (1.46) 
whereas the others let us write (jl.41j) as 

hjk = (—) tu 2 D jk e^ r -^ (1.47) 



h TT = 





h TT = 

11 'XX 


-h TT 


h TT = 

11 'xy 


V r 



We have still freedom of adopting TT gauge and may use it to achieve further 
simplification. Let us choose z axis along the direction of propagation of the 
wave. We will then have 

(1.48) 

= -u 2 {l xx -l yy ) — (1.49) 
u% y e iujr (1.50) 

where Ijk represent the trace free part of quadrupole moment tensor i.e. 

l jk = I jk - U ]k l\ (1.51) 

As an application of our results let us determine the amplitude of GW gener- 
ated by a laboratory source. 

Laboratory generator (Bar) 

Let us take a system of two equal mass points capable of oscillating about 
their mean position, such as shown in Figure (|1.3j) . Let the system oscillates 
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longitudinally with angular frequency w about its mean position i.e. 

X\ = — \l — Acoswt ] 

1 2 \ (1.52) 

X2 = ^/o + Acoswt j 

where l a is the normal separation between the mass points and A represents 

the amplitude of oscillation. Now it is straight-forward to compute Ijk- The 

only non-zero component is 

I xx = m[(x 1 ) 2 + (x 2 ) 2 ] 

= const. + mA 2 cos 2wt + 2m/ D A cos wt 

For purpose of wave generation the constant term is irrelevant. We may obtain 

the non- vanishing components of Ij^ as 



1.53) 



I xx = |m/ Ae- iw * + |m/ A 2 e- 2iwt 

I yy = \ zz = -| m / Ae- iw * - |m/ A 2 e- 2iwt 



1.54) 



hll = -Ky = -[2mw 2 / Acos(w(r-t)) 



Finally one obtains after taking the real part 

- [2mw 2 / D A 

4mw 2 A 2 cos(2w(r-t))]/r } (1-55) 
For a laboratory generator, let us take 



^ = 



m=10 6 kg, l = lm, A = 10 m , w = 10" 4 s' 1 (1.56) 

The data chosen represents a typical bar detector. Substituting the values 
after converting them into geometrised units (G = 1 = c) the amplitude of 
generated wave is about 10~ 34 /r 

\h\ ~ 10~ 34 /r ; Laboratory source (1-57) 

This demonstrates that laboratory generators are unlikely to produce use- 
ful GW for its demonstration. For sake of comparison, let us conclude this 
subsection by making an estimate of the strength of waves from strongest 
astrophysical sources. 
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Figure 1.3: An idealised gravitational wave detector with two identical masses 
connected with a spring. 

Astrophysical sources 

For a strong GW we should have h^ u = 0(1). This would occur near the source 
where the Newtonian potential would be of the order 1. For a source of mass 
M this should be at a distance of order M. As we have seen the amplitude of 
GW falls off as r _1 far from the source. This means that the largest amplitude 
expected to be incident on earth would be ~ M/R where R is the distance of 
the source from earth. For the formation of 10 M Q black hole in a supernova 
explosion in a nearby galaxy 10 23 away this is about 1CT 17 . Thus 

|/i| maa , ~ 10~ 17 ; Astrophysical sources (1.58) 

Thus is in fact an upper limit and less violent events will lead to very much 
smaller amplitudes. 



1.2 Gravitational wave detectors 

Equation (jl.HOj) may be recast as 

= (1-59) 

This means that the strain produced due to the hitting of a GW in two mass 
points separated by e is of the order of h. In view of Eq. fjl.58|) let us note that 
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even a strong GW signal would produce a space strain of the order of 1CT 17 , 
an unbelievably small effect which would jerk masses spaced at one kilometer 
by a mere 10~ 20 -one thousand of the diameter of a proton! 

Josheph Weber (1960) who pioneered the direct detection of GW con- 
structed an instrument consisting of massive cylinders of aluminium so-called 
"bar" detectors. These detectors exploit the sharp resonance of the cylinder 
to get their sensitivity which is normally confined to a narrow bandwidth (one 
or a few Hz) around the resonant frequency. 

Despite their great potential sensitivity the primary drawback of resonant 
bars is that they are by definition resonant. They are sensitive mainly to the 
signal with a frequency corresponding to the bar mechanical ringing frequency 
of the order of 1 kHz. A bar would respond to the hammer blow of the 
asymmetrical supernova explosion by simply ringing at its own bell tones and 
would be excited by a twin neutron star inspiral only in that brief instant 
when the two stars crib up through the bell tone frequency. 

Bar detectors continue to be developed, and have until very recently had a 
sensitivity to broadband bursts. However, the best hope for the first detection 
of GW lies with large-scale interferometers. Within ten years from now, we 
may see the launch of space-based interferometer, LISA, to search for signals at 
frequencies lower than those that are accessible from the ground. To measure 
the strain produced by GW to bar or interferometric detector, one must fight 
against the different sources of noises. 

An interesting additional twist is given by the fact that the gravitational 
wave may be accompanied by gamma ray bursts. The GW detectors will then 
work in coincidence not only with themselves and GW bar antennas, but also 
with conventional high-energy physics detectors like the underground neutrino 
experiments and the orbital gamma rays burst monitors. 
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In this section we review, in brief, these developments. In subsection (jl.2.11) 
the underlying physical principles of the bar detectors have been described. 
Subsection (jl.2.2j) deals with large-scale interferometers. Finally subsection 
(jl.2.3|) describe the space-based project LISA, the exciting prospects for the 
long-term future of GW detection in space. For details, one may refer to 
Thorne (1987, 1995), Saulson (1994). 

1.2.1 Bar detectors 

A bar detector, in its simplest form, may be idealised as a system of two mass 
points coupled to a spring and may be depicted as shown in Fig. (jl.3|) . Let 
the system lies on the x axis of our TT coordinate system with the masses at 
coordinate positions X\ and x 2 . The force free oscillation of the system, in flat 
space time, may be expressed via 

mii,oo = -k(x! - x 2 + l ) - v{x\ -x 2 ) jG 
ma^oo = -k(x 2 - x x + l ) - v{x 2 - Zi) , 

where l Q , k and v represent respectively unstretched length of the spring, spring 
constant and damping constant. We can combine these equations to obtain 
the usual damped harmonic oscillator equation 

Coo + 2 7 f, + «& = (1.61) 

by introducing 

£ = x 2 — x\ — l Q , wl = 2n/m, j — f/m (1-62) 

Let us recall that the TT coordinate system is not convenient for discussion 
of the dynamics of such system because in this frame a free particle always 
(before the arrival and after the passage of the wave) remains at rest. However, 
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this fact is useful in assigning a local inertial frame {x a } at some TT coor- 
dinate. Suppose that the only motions in the system are those produced by 
the wave then masses velocities will be very small and we may apply Newton's 
equations of motion for the masses: 



mx J , , = F j 



1.63) 



where {F J '} are the components of any non gravitational forces on the mass. 
Further as the coordinates {a;- 7 } differ negligibly to the order of from that 
of its value in TT coordinate system, we may write this equation with 
negligible error: 



mar 



oo 



1.64) 



The only non gravitational force on each mass is due to the spring. The 
spring will exert a force proportional to its instantaneous proper extensions. 
If the proper length of the spring is I and the direction of propagation of the 
wave, for simplicity, is assumed to be along z axis then 



l(t) 



r x 2(t) 

/ [1 + /&"(*)] 1/2 d* 



i + (*) 



lx 2 - X\, 



1.65) 



(refer to (|1.30p ). Hence, equation of motion of the system after hitting of the 



wave is given via 



ma^oo = -k(1 -1) - v(l — l) ,o 
mx 2 ,oo = -l ) - v{l -l ) ,o 



Let us define 



leading to 



x 2 - Xi 



i + if I 



(x 2 - Xi) - l 



1.66) 



1.67) 



{i + Q (l - \h T x f\ =i + l - \h T Zl + O (| VI 2 ) (1.68) 
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Using this equation, we may obtain from Eq. (jl.66|) 

Coo + 2 7 £,o + ^=^oCoo (1-69) 

correct to first order in h^J. This is the fundamental equation governing the 
response of the detector to the GW. It has the simple form of a forced, damped 
harmonic oscillator . 

Let a GW of frequency u described via 

h^J = A cos tut (1.70) 

hits the detector then the steady solution for £ may be taken as 
£ = RcosfW + <f) ; 

R = \1 u 2 K/[{w -uj) 2 + ^ 2 v 2 } 1/2 , } (1.71) 
tanip = 2ulu/ (w 2 — uo 2 ) 
The average energy of oscillation of the detector over one period, 2tt/uj: 

(E) = ^mR 2 (w 2 + u 2 ) (1.72) 

If we wish to detect a specific source whose frequency u is known, then we 
should adjust w Q equal to u for maximum response (resonance). The resonance 
amplitude and energy of the detector will be 

^■resonant = ^(A-^/l (1-73) 

x 

^resonant = — mlfe 2 A 2 (u / (1.74) 

64 

The ratio (w/7) is related to what is called the quality factor Q. 

Q = cu/2 7 (1.75) 

Kesonant = -^Hl/ 2 ^ 2 A 2 Q 2 (1-76) 

lb 

The bar detectors are massive cylindrical bars; its elasticity provides the 
function of the spring. When waves hit the bar broadside, they excite its longi- 
tudinal modes of vibration. The first detectors built by Weber were aluminium 



20 



bars of mass 1.4 x 10 3 kg, length l = 1.5 m, resonant frequency w Q = 10 4 s _1 
and Q about 10 5 . This means that a strong resonant GW of A = 10~ 20 will 
excite the bar to an energy of the order of 10~ 20 J. The resonant amplitude 
given by (jl.73J) is only about 10~ 15 m, roughly the diameter of an atomic 
nucleus. 




Figure 1.4: A schematic of resonant bar detector. 

Clearly, the detections of such small levels of excitations will be hampered 
by random noise in the oscillator. For example, thermal noise in any oscillator 
induces random vibration with a mean energy of kT; k the Boltzmann constant 
having value 1.38 x 10~ 23 J/K. At room temperature (T ~ 300K) the thermal 
noise amounts to the energy ~ 4 x 10~ 21 J. Other sources of noise such as 
vibrations from passing vehicles and every day seismic disturbances could be 
considerably larger than this, so the apparatus has to be carefully isolated. 

A typical "bar" detector consists of a cylinder of aluminium with a length 
l Q ~ 3 m, a resonant frequency of order w Q ~ 500 Hz to 1.5 kHz, and a mass 
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Detector 


Location 


Taking data since 


NAUTILUS 


Frascati, Rome 


1993 


EXPLORER 


Cern (Rome group) 


1990 


ALLEGRO 


Louisiana, USA 


1991 


AURIGA 


Padua, Italy 


1997 


NIOBE 


Perth, Australia 


1993 



Table 1.1: The existing resonant bar detectors. 

~ 1000 kg whose mechanical oscillations are driven by GW [Fig. (jl.4j) ]. A 
transducer converts mechanical vibrations of the bar into an electrical signal, 
which is then amplified by an amplifier and recorded. 

Currently there are a number of bar detectors in operation [see Table fll.ljl ]. 
Some of these operate at room temperature and some others at cryogenic tem- 
perature. Some detectors (NAUTILUS & EXPLORER) may be cooled down 
to ultra cryogenic temperature. They can detect signal amplitudes h ~ 10 -20 
in a band width of 10-20 Hz around a central frequency of 1 kHz. Asymmetric 
supernovae in our Galaxy are the best candidates for these detectors. For 
example, a supernova collapse in our galaxy at a distance of 10 kpc emits GW 
of the amplitude h ~ 10~ 17 . At present this sensitivity has been achieved by 
some of the bar detector. They may also see continuous radiation emitted by 
a neutron star if the frequency happens to lie in their sensitivity band. 
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1.2.2 Ground based laser interferometric detectors 

The effect of GW is to produce a transverse shear strain and this fact makes 
the Michelson interferometer an obvious candidate for a detector. The Michel- 
son interferometers must have kilometric arms, constituted by "high fineness" 
Fabry Perot cavities to trap the light for long period and to increase the sensi- 
tivity. Laser standing power measured in KW will be stored within the cavities. 
Beam losses at the level of 10~ 6 per passages are required. Mirrors must be 20 
or 30 cm in diameter just to hold diffraction losses, and coating with parts per 
million reflection losses have been developed. Vacuum pipe with diameters of 
the order of a meter are necessary just to contain the diffraction limited laser 
beams and vacuum exceeding 10 torr is required to prevent the introduction 
of spurious signal by refraction index fluctuations. The list of ground based 
laser interferometric detector site and orientation are shown in Table (jl.2j) . 

Interferometric detectors currently under construction will increase our 
ability to directly observe GW. The Japanese TAMA-300 meter interferometric 
detector having sensitivity 10~ 20 is already operational and some initial hours 
of data has been analysed. Several other projects are now nearing to comple- 
tion: The British-German collaboration is constructing a 600 m interferometer 
(GEO) in Hannover, Germany, the French-Italian collaboration is building a 3 
km detector (VIRGO) near Pisa, Italy and the Americans are building two 4 
km antennas (LIGO), one in Livingston and the other in Hanford in the USA 
which will have the sensitivity 10~ 21 . These detectors will start taking data 
between 2002 and 2003. The larger of these detectors, LIGO and VIRGO, are 
likely to be upgraded in sensitivity by an order of magnitude with a better 
low-frequency performance in 2005. These ground based interferometers will 
eventually be sensitive to sources in the frequency range from 10 Hz to several 
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Detector 


Location 


Length (m) 


Corner Location 


Arm 1 


Arm 2 


Glasgow 


Glasgow, GBR 


10 


55.87°N - 4.28°W 


77.00° 


167.00° 


CIT 


Pasadena, CA, USA 


40 


34.17°N - 118.13°W 


180.00° 


270.00° 


MPQ 


Garching, GER 


30 


48.24°N - 11.68°W 


329.00° 


239.00° 


ISAS-100 


Tokyo, JPN 


100 


35.57°N - 139.47°W 


42.00° 


135.00° 


TAMA-20 


Tokyo, JPN 


20 


35.68°N - 139.54°W 


45.00° 


315.00° 


Glasgow 


Glasgow, GBR 


10 


55.87°N - 4.28°W 


62.00° 


152.00° 


TAMA-300 


Tokyo, JPN 


300 


35.68°N - 139.54°W 


90.00° 


180.00° 


GEO-600 


Hannover, GER 


600 


52.25°N 9.81°W 


25.94° 


291.61° 


VIRGO 


Pisa, ITA 


3000 


43.63°N 10.5°W 


71.50° 


341.50° 


LIGO 


Hanford, WA, USA 


4000 


46.45°N 119.41°W 


36.80° 


126.80° 


LIGO 


Livingston, LA, USA 


4000 


30.56°N 90.77°W 


108.00° 


198.00° 



Table 1.2: Site and orientation of Earth based interferometric gravitational 
wave detectors. 

kHz. 

When a GW, with frequencies high compared to the masses ~ 1 Hz pen- 
dulum frequency, passes through the detector of arm length l Q then one arm 
of the detector get stretched in one direction whereas the other arm gets com- 
pressed. If the total change in the the arm-length is 51, then the response R(t) 
of the detector which is defined as f- may be given via 



R(t) = T = F + h + (t) + F x h x (t) (1.77) 

where, h + , h x are two polarisations of the signal and F + , F x are beam pattern 
functions arising because of transformation equations relating to the wave 
frame and the detector frame [refer to Eq. (jl.31j) ]. Let the detector frame 
(x, y, z) and the wave frame (X, Y, Z) are related via transformation matrix 
A i.e 

e = A* R e, e = A k ^ (i.78) 

where capital latin letters refers to wave frame quantities and lower case for the 
detector frame quantities. It can be shown that the response of the detector 
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may be expressed as 
A/ 

- = - sin 2ft [(A x x A y Y - A X Y A Y ) h+ + (A X X A V Y + A Y A y x ) h x ] (1.79) 

'o 



where 2ft is the angle between the arms of the detector, normally kept equal 
to 90° [Schutz and Tinto (1987)]. 




Figure 1.5: A schematic of Earth based laser interferometric gravitational wave 
detector. 

A laser interferometric detector consists of four test masses/mirrors that 
hang from vibration isolated supports as shown in Fig. (jl.5|) . The mirrors 
facing each other on each arm form a Fabry-Perot cavity. A beam splitter 
splits a carefully prepared laser beam in two, and directs the resulting beams 
down the two arms. Each beams penetrates through the mass, and through 
the dielectric coating (the mirror). The length of the arm's of the Fabry-Perot 
cavity is adjusted to be nearly an integral number of half wave-length of light. 
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The cavity's end mirror has much higher reflectivity than its corner mirror, so 
the trapped light leaks back out through the corner mirror, and then hits the 
beam splitter where it recombines with light from the other arm. Most of the 
recombined light goes back toward the laser (where it can be returned to the 
interferometer by a "power recycled mirror" labeled PRM, but a tiny portion 
goes toward the photo-diode. 

When a GW hits the detector and moves the masses, thereby changing the 
lengths of the arms of the two cavities, it shifts each cavity's resonant frequency 
slightly relative to the laser frequency, and the phase of the light that exists 
from the cavity toward the beam splitter. Correspondingly, the relative phase 
of the two beams returning to the splitter is altered by an amount A$ oc 51, 
and this relative phase shift causes a change in the intensity of the recombined 
light at photo-diode, AI pd oc A$ oc 51 oc R(t). Thus, the change of photo- 
diode output current is directly proportional to the GW strain. 

Scientists are now studying the technologies that may be needed for a fur- 
ther improvement of the sensitivity of the next-generation detectors. This may 
involve cooling mirrors, using ultra-massive mirrors of special materials, us- 
ing purely non-transmissive optics, and even circumventing the quantum limit 
in interferometers, as has been studied for bars. The goal of next-generation 
detectors would be limited just by gravity-gradient noise and quantum effects. 

1.2.3 Laser interferometric space antenna 

Gravity-gradient noise on Earth is much larger than the amplitude of any 
expected waves from astronomical sources at frequencies below about 1 Hz, 
but this noise falls off a 1/r 3 as one moves away from the Earth. A detector 
in space would not notice the Earth's noisy environment. The Laser inter- 
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ferometric space antenna (LISA) project, efforts underway by both European 
Space Agency (ESA) and NASA with a view toward a collaborative mission 
around 2010, would open up the frequency window between 0.1 mHz and 0.1 
Hz. 

5 x 10 — 




Figure 1.6: A schematic diagram of the laser interferometric space antenna. 

A concept of the project is shown in Fig. (jl.6|) . Three spacecraft are placed 
in solar orbit at 1 A.U, about 20 degrees behind the Earth in it orbit. The 
spacecraft are located at the corners of an equilateral triangle with 5 x 10 6 
km-long sides. Two arms of the triangle comprise a Michelson interferometer 
with vertices at the corners. The third arm permits another interferometric 
observable to be measured, which can determine a second polarization. The 
interferometers use one micron light as the terrestrial detectors but need only 
a single pass in the arms to gain the desired sensitivity. The end points of 
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the interferometers are referenced to proof masses free-floating within and 
shielded by the spacecraft. The spacecraft is incorporated in a feedback loop 
with precision thrust control to follow the proof masses. 

The main environmental disturbances to LISA are the forces from the Sun.: 
fluctuations in solar radiation pressure and pressure from the solar wind. To 
minimize these, LISA incorporates drag-free technology. Interferometry is ref- 
erenced to an internal proof mass that falls freely, not attached to the space- 
craft. The job of the spacecraft is to shield this mass from external distur- 
bances. It does this by sensing the position of the mass and firing its own jets 
to keep itself (the space craft) stationary relative to the proof mass. To do this, 
it needs thrusters of very small thrust that have accurate control. The avail- 
ability of such thrusters, of the accelerometers needed to sense disturbances 
to the spacecraft, and of lasers capable of continuously emitting 1 W infrared 
light for years, have enabled the LISA mission. 

LISA is supposed to see many exciting sources for example the coalescences 
of giant black holes in the centre of galaxies. LISA will see such events with ex- 
traordinary sensitivity, recording typical signal-to-noise-ratios of 1000 or more 
for events at redshift 1. 

1.3 Sources of gravitational waves 

Astronomical observations have led to the belief that luminous matter consti- 
tutes a small fraction of the total matter content of the Universe. More than 
90% of the mass in the Universe is electro-magnetically silent. The presence 
of dark matter is inferred from the gravitational influence it causes on lumi- 
nous matter. It is possible that some fraction of this dark matter is a strong 
emitter of GW. There are many reviews on GW sources (Thorne, 1987; Blair, 
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1991; Schutz, 1989, 1993, 1999, Sathyaprakash, 1999 and Grishchuk, et. al., 
2000). The discussion on the GW sources in this section is introductory and 
for details one may refer to these reviews. The anticipated GW sources can 
be classified into (i) Transients (ii) Continuous and (iii) Stochastic. 

1.3.1 Transient sources 

Supernovae and inspiralling compact stars in binary systems are prime ex- 
amples of sources of this type. These sources last for a duration of about 
10~ 3 - 10 4 seconds. 

Supernovae explosions 

The type II Supernovae explosions, which are believed to occur as a result of 
the core collapse of an evolved massive (> 9M ) star and which are associated 
with violent mass ejection with velocities of order 0.03 with formation of a 
compact remnant — neutron star or a black hole-may emit significant amount 
of GW depending upon how asymmetric the collapse is. The emission oc- 
curs essentially during the rotational core collapse, bounce and oscillations, 
rotation-induced bars and convective instabilities set up in the core of the new 
born neutron star. 

Rapid rotation flattens the collapsing core inducing large quadrupole mo- 
ment; thus generating GW. Study of a wide range of rotational core collapse 
models suggests that the largest signals are produced by models which are 
(i) initially slowly rotating and have a stiff equation of state or (ii) initially 
rapidly rotating and have a soft equation of state. In the first case bounce 
occurs at densities above nuclear matter density, with a fast deceleration of 
the collapsing core resulting in the emission of GW signals. In the second 
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case, the quadrupole moment is large due to rapid rotation which facilitates 
emission of GW. However, in either case the signals are not strong enough to 
be interesting sources for the first generation detectors. 

When the core's rotation is high enough it may cause the core to flatten 
before it reaches nuclear density leading to an instability which transforms the 
flattened core into a bar-like configuration spinning about its transverse axis. 
Some of these instabilities could also fragment the core into two or more pieces 
which then rotate about each other. Both are efficient ways of losing energy 
in the form of GW. 

Instabilities in the core of the newly born neutron star, which last for about 
a second after the collapse, are likely to produce GW due to anisotropic mass 
distribution and motion. 

Inspiraling compact binaries 

The binary systems whose either member is a compact star e.g. neutron star 
(NS) and black hole (BH) are the most promising transient sources of GW 
during the phase of their coalescence. The well known binary pulsar PSR 
1913 + 16 is such a system but it will coalesce in a time scale of 10 9 years from 
now — not a right candidate. However, there are binaries in our galaxy with 
coalescence time scale much shorter than this. Further, statistical analysis of 
binary pulsars estimates three NS-NS coalescence per year out to a distance of 
200 Mpc (Phinney, 1991; Narayanan et. al 1991). The initial LIGO/VIRGO 
interferometers have a fair chances to see the inspiral events. Compact inspi- 
raling binaries emit quasi-periodic GW with a frequency that sweeps upward 
toward a maximum frequency. The maximum frequency may be of the order 
of 1 KHz for neutron stars. In the lower frequency regime the wave form is 
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easily computed from the quadrupole formalism. At higher frequencies Post- 
Newtonian corrections would be required (Thorne, 1987; Krolak, 1989). In 
view of the strong potentialities of such sources the various aspects related to 
the emission of G W has been dealt extensively in full details and one may refer 
to Sathyprakash, 1999 and Grischuk, et al. 2000. 

1.3.2 Continuous gravitational wave 

The continuous gravitational wave (CGW) sources are of prime interest be- 
cause such sources can be observed again and again and hence single inter- 
ferometer is sufficient to confirm its detection. However, one can't expect a 
source to be continuous as well as strong. For emission of GW from pulsars, 
there should be some asymmetry in it. There are several mechanisms which 
may lead to deformations of the star, or to precession of its rotation axis. The 
characteristic amplitude of GW from pulsars scales as 




where I is the moment of inertia of the pulsar, f is the GW frequency, e is a 
measure of deviation from axissymmetry and r is the distance of the pulsar. 

As remarked earlier, pulsars are born in supernovae explosions. The outer 
layers of the star crystallizes as the newborn star pulsar cools by neutrino 
emission. Estimates, based on the expected breaking strain of the crystal 
lattice, suggest that anisotropic stresses, which build up as the pulsar looses 
rotational energy, could lead to e < 10~ 5 ; the exact value depends on the 
breaking strain of the neutron star crust as well as the neutron star's "geolog- 
ical history" , and could be several orders of magnitude smaller. Nonetheless, 
this upper limit makes pulsars a potentially interesting source for kilometer 
scale interferometers. 



31 



Large magnetic fields trapped inside the super fluid interior of a pulsar 
may also induce deformation of the star. This mechanism has been explored 
recently, indicating that the effect is extremely small for standard neutron star 
models ( e < 1CT 9 ). 

Another plausible mechanism for the CGW is the Chandrasekhar-Friedman- 
Schutz (CFS) instability, which is driven by GW back reaction. It is possible 
that newly-formed neutron stars may undergo this instability spontaneously 
as they cool soon after formation. The frequency of the emitted wave is deter- 
mined by the frequency of the unstable normal mode, which may be less than 
the spin frequency. 

Accretion is another way to excite neutron stars. There is also the Zim- 
mermann - Szedinits mechanism where the principal axes of moment inertia 
are driven away from the rotational axes by accretion from a companion star. 
Accretion can in principle produce relatively strong wave since the amplitude 
is related to the accretion rate rather than to structural effects in the star. 

1.3.3 Stochastic waves 

Catastrophic processes in the early history of Universe, as well as the as- 
trophysical sources distributed all over the cosmos, generate stochastic back- 
ground of GW. A given stochastic background will be characteristic of the 
sources that are responsible for it and it may be possible to discriminate dif- 
ferent backgrounds provided their spectral characteristics are sufficiently dif- 
ferent . 

Primordial background: It is believed that similar to cosmic microwave 
background (CMBR) the GW background was also produced at the same 
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time as a result of quantum fluctuations in the early Universe (Grischuk, 1997). 
Primordial background radiation would freely travel to us from almost the very 
moment of creation because GW couples very weakly with matter. Hence, its 
detection would help us to get a picture of the first moments after the big 
bang. The COBE data have set limits on the strength of the GW background. 
The strength are far too weak to be detected by any planned ground-based 
detectors. However, advanced LIGO detectors may observe the background 
generated by the collisions of a cosmic string network (for details one may 
refer to Allen, 1997). 

Phase transitions in the early Universe, inspired by fundamental particle 
physics theories, and cosmological strings and domain walls are also sources 
of a stochastic background. These processes are expected to generate a back- 
ground which has a different spectrum and strength than the primordial one. 
Future ground-based detectors will achieve good enough sensitivity to measure 
this background GW and such measurements will prove to be a good test bed 
for these cosmological models. 

Supernovae background: Even though an individual supernova may not be 
detectable out to a great distance, the background produced by all supernovae 
within the Hubble radius might indeed be detected (Blair & Ju, 1996). The 
coincidence between an advanced interferometer and a resonant bar within 50 
km of the interferometer will enable the detection of this background. These 
studies may shed light on the history of star formation rate, a subject of vig- 
orous debate amongst astrophysicists. 

Galactic binary background: Binaries with orbital periods P ~ 10~ 4 — 10~ 2 s 
will be observable in space-based detectors. A large number of them are present 
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in our galaxy but they will not be identifiable separately because they are at a 
large distances and consequently have a feeble amplitude. However they would 
contribute to background GW. In addition to binaries of compact stars, there 
are also other binaries consisting of white dwarfs, cataclysmic variables, etc., 
which will also contribute to the background radiation produced by compact 
binaries. The net effect is that these sources will appear as a background noise 
in space interferometers. By studying the nature of this background one can 
learn a lot about binary population in our Galaxy. 

Galactic pulsar background The CGW from pulsars could also produce 
background of radiation which will limit the sensitivity of the ground-based 
laser interferometers (Giazotto, Gourgoulhon and Bonazzola, 1997). There are 
about 10 9 neutron stars in our Galaxy of which about 2 x 10 5 will contribute to 
the background radiation. This background radiation will be prominent, and 
observable in the LIGO/VIRGO detectors, in the frequency range of 5-10 Hz, 
at an rms amplitude of h rms ~2x 10~ 26 , where the rms is computed over 10 5 
sources. At frequency of 10 Hz, wavelength of GW will be around 30,000 Km. 
Hence, it would be possible to cross-correlate data from two distant detectors, 
such as two LIGOs, or LIGO and VIRGO, and discriminate the background 
against other sources of noise. 
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Chapter 2 

Data analysis concept 

2.1 Introduction 

The GW data analysis strategy is different in many ways from conventional 
astronomical data analysis. This is due to the following: 

• GW antennas are essentially omni-directional with their response better 
than 50% of the average over 75% of the sky. Hence our data analysis 
systems will have to carry out all-sky searches for the sources. 

• Interferometers are typically broad-band covering 3 to 4 orders of mag- 
nitude in frequency. While this is obviously to our advantage, as it helps 
to track sources whose frequency may change rapidly, it calls for searches 
to be carried over a wide-band of frequencies. 

• In contrast to electromagnetic (EM) radiation, most astrophysical GWs 
are tracked in phase and the signal-to-noise ratio (SNR) is built up by 
coherent superposition of many wave cycles emitted by a source. Conse- 
quently, the SNR is proportional to the amplitude and only falls off, with 
the distance to the source r, as 1/r. Therefore, the number of sources of 
a limiting SNR increases as r 3 for a homogeneous distribution of sources 
in a flat Universe, as opposed to EM sources that increase only as r 3 / 2 . 
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• The GW antennas acquire data continuously for many years at the rate of 
several mega-bytes per second. It is expected that about a hundredth of 
this data will have to pass through our search analysis systems. Unless 
on-line processing can be done we cannot hope to make our searches. 
This places huge demands on the speed of our data analysis hardware. 
A careful study of our search algorithms with a view to making them 
as optimal (maximum SNR) and efficient (least search times) as one 
possibly can is required. 

2.2 Gravitational wave antenna sensitivity 

The performance of a GW detector is characterised by the one sided power 
spectral density (PSD) of its noise background. The analytical fits to noise 
power spectral densities S n (f) of ground based interferometers are given in 
Table (|2.1|) where S Q and f no represent respectively the value of minimum noise 
and the corresponding frequency. At the lower-frequency cutoff f] and the 
high-frequency cutoff f u , S n (f) can be treated as infinite. One can construct 
the noise PSD as follows: 

A GW detector output represents a dimensionless data train, say x{t). In 
the absence of any GW signal. The detector output is just an instance of noise 
n(t), that is, x(t) = n(t). The noise auto-correlation function c is defined as 



where ( ) represents the average over an ensemble of noise realisations. In 
general, c depends both on ti and t 2 - However, if the detector output is a 
stationary noise process, i.e. its performance is, statistically speaking, inde- 
pendent of time, then c depends only on r = t 2 — We shall, furthermore, 
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Detector 


Fit to noise power spectral density 


S 
(Hz)- 1 


fn 

(Hz) 


/i 
(Hz) 


fu 

(Hz) 


LIGO I 


Sol* \{fn /f) i + 2(f/f n ) 2 1 


8.0 x 10~ 46 


175 


40 


1300 


LIGO II 


5 /ll {2(/„ //) 9 / 2 + 9/2[l + (///„J 2 ]} 


7.9 x 10~ 48 


110 


25 


900 


LIGO III 


So/5 {(/„ o //) 4 + 2[l + (///n ) 2 ]} 


2.3 x 1(T 48 


75 


12 


625 


VIRGO 


S / A [290(/ s //) 5 + 2(/„ o //) + 1 + (///„J 2 ] 


1.1 x 10~ 45 


475 


16 


2750 


GEO600 


W5 [4(/„ o //)3/ 2 -2 + 3(///„J 2 ] 


6.6 x 1(T 45 


210 


40 


1450 


TAMA 


So/32 {{f n Jff + 13(/ no //) + 9[1 + (f/.f no ) 2 }} 


2.4 x 10~ 44 


400 


75 


3400 



Table 2.1: Analytical fits to noise power spectral densities S n (f) of ground 
based interferometers. 
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assume that c(r) = c(— r). For data from real detectors the above average can 
be replaced by a time average under the assumption of ergodicity: 

1 f T ' 2 

c(r) = -/ n(t)n(t-r)dt. (2.2) 

The assumption of stationarity is not strictly valid in the case of real GW 
detectors; however, if their performance does not vary greatly over time scales 
much larger than typical observation time scales, stationarity could be used 
as a working rule. While this may be good enough in the case of binary 
inspiral and coalescence searches, it is a matter of concern for the observation 
of continuous and stochastic GW. Under such an assumption the one-sided 
noise PSD, defined only at positive frequencies, is the Fourier transform (FT) 
of the noise auto-correlation function: 

Sn{f) = 2 / C ( r ) e ~ 27 " /Trfr > / ^ > 

= 0, /<0, (2.3) 

where a factor of 1/2 is included by convention because it has been assumed 
that c(r) is an even function. This equation implies that S n (f) is real. It is 
straightforward to show that 

(n(f)n m (f)) = S n (f)S(f-f), (2-4) 

where n(f) represents the Fourier transform of n(t) and * denotes complex 
conjugation. The above identity implies that S n (f) is positive definite. 

The autocorrelation function c(r) at r = can be expressed as an integral 
over S n (f). Indeed, it is easy to see that 

POO 

(n 2 (t)) = 2 / S n (f)df. (2.5) 
Jo 

The above equation justifies the name power spectral density given to S n (f). 
It is obvious that S n (f) has dimensions of time but it is conventional to use, 
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instead, the dimensions of Hz -1 since it is a quantity denned in the frequency 
domain. The square-root of S n (f) is the noise amplitude, \J S n (f), and has 
dimensions of per root Hz. It is often useful to define the dimensionless quantity 
hn(f) = fS n (f), called the effective noise. In GW interferometer literature 
one also comes across the displacement noise or strain noise defined as hi a (f) = 
l h n (f), and the corresponding noise spectrum Si a (f) = l 2 S n (f), where l is 
the arm length of the interferometer. The displacement noise gives the smallest 
strain Sl/l Q in the arms of an interferometer which can be measured at a given 
frequency. 

Sensitivity vs source amplitudes 

One compares the GW amplitudes of astronomical sources with the instru- 
mental sensitivity and assesses what sort of sources will be observable in 
the following way. Firstly, as comparisons are almost always made in the 
frequency-domain it is important to note that the Fourier component h(f) of 
a deterministic signal h(t) has dimensions of Hz -1 and the quantity f\h(f)\, 
is dimensionless. It is this last quantity that should be compared with h n (f) 
to deduce the strength of a source relative to detector noise. Secondly, it is 
quite common also to compare the amplitude spectrum per logarithmic bin 
of a source, VJ\h(f)\, with the amplitude spectrum of noise, a/ S n (f), both 
of which have dimensions of per root Hz. For monochromatic sources, one 
compares the effective noise in a long integration period with the expected 
"instantaneous" amplitudes in the following way: A monotonic wave of fre- 
quency f observed for a time T Q is simply a narrow line in a frequency bin 
of width Af = 1/T Q around f . The noise in this bin is S n (f)Af = S n {f)/T . 
Thus the SNR after a period of observation T Q is 

— = ^° (2 6) 
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One, therefore, computes this dimensionless noise spectrum for a given du- 
ration of observation, S n (f)/T , to assess the detectability of a continuous 
GW. Surely, if the observation time is T Q then the total energy (that is, the 
integrated power spectrum) of both the signal and noise must increase in pro- 
portion to Tg. Then how does the SNR for a continuous wave improve with the 
duration of observation? The point is that while the signal energy is all con- 
centrated in one bin, the noise is distributed over the entire frequency band. 
As T Q increases, the frequency resolution improves as 1/T G and the number of 
frequency bins increase in proportion to T Q . Consequently, the noise intensity 
per frequency bin increases only as T a . Now, the signal intensity is concentrated 
in just one bin since the signal is assumed to be monochromatic. Therefore, 
the power SNR increases as T G , or, the amplitude SNR increases as \f 7 T . 

2.3 Noises in the earth based interferometric 
detectors 

As mentioned in the earlier section the performance of GW detector is char- 
acterised by the one-sided PSD. The sensitivity of ground based detectors is 
limited at frequencies less than a Hertz by the time- varying local gravitational 
field caused by a variety of different noise sources, e.g. low frequency seismic 
vibrations, density variation in the atmosphere due to winds, etc. Thus, for 
data analysis purposes, the noise PSD is assumed to be essentially infinite 
below a certain lower cutoff f s . The effective noise h n (f) expected in these 
detectors is equal to \/fS n (f). Type of noises that one have to fight to detect 
the signals are: 
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Ground vibration: External mechanical vibrations are a problem for bar 
detectors, too, but are more serious for interferometers, not least because in- 
terferometers bounce light back and forth between the mirrors, and so each 
reflection introduces further vibrational noise. Suspension/isolation systems 
are based on pendula. A pendulum is a good mechanical filter for frequencies 
above its natural frequency. By hanging the mirrors on pendula of perhaps 
0.5 m length, one achieves filtering below a few Hz. Since the spectrum of 
ground noise falls at higher frequencies, this provides suitable isolation. But 
these systems can be very sophisticated; the GEO600 detector has a three- 
stage pendulum and other vibration isolation components. The most ambi- 
tious isolation system is being developed for the VIRGO detector. 

Thermal noise: Vibrations of the mirrors and of the suspending pendu- 
lum can mask gravitational waves. As with vibrational noise, this is increased 
by the bouncing of the light between the mirrors. Unlike bars, interferometers 
measure only at frequencies far from the resonant frequency, where the am- 
plitude of vibration is smaller. Thus, the pendulum suspensions have thermal 
noise at a few Hz, but measurements will be made above 20 or 30 Hz in the first 
detectors. Internal vibrations of the mirrors have natural frequencies of several 
kilohertz. By ensuring that both kinds of oscillations have very high Q, one can 
confine most of the vibration energy to a small bandwidth around the resonant 
frequency, so that at the measurement frequencies the vibration amplitudes 
are small. This allows interferometers to operate at room temperature. But 
mechanical Qs of 10 7 or higher are required, and this is technically demanding. 

Shot noise: The photons that are used to do interferometry are quan- 
tized, and so they arrive at random and make random fluctuations in the light 
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intensity that can look like a gravitational wave signal. The more photons 
one uses, the smoother will be the interference signal. As a random process, 
the error improves with the square-root of the number N of photons. Using 
infrared light with a wavelength A ~ 1 /an, one can expect to measure to an 
accuracy of 

Ah shot ~ A/(2ttv / N) 

To measure at a frequency /, one has to make at least 2f measurements per 
second i.e. one can accumulate photons for a time 1/2/. With light power 
P and GW amplitude h one gets N = P/(hc/X)/(2f) photons. In order that 
AL S hot should be below 1CT 16 m one needs large light power, far beyond the 
output of any continuous laser. 

Light-recycling techniques overcome this problem by using light efficiently 
An interferometer actually has two places where light leaves. One is where 
the interference is measured. The other goes back towards the input laser. 
Normally one arranges that no light goes to the interference sensor, so that 
only when a gravitational wave passes does a signal register there. This means 
that all the light normally returns to the mirror, apart from small losses at the 
mirrors. Since mirrors are of good quality, only one part in 10 3 or less of the 
light is lost during a 1 ms storage time. By placing a power-recycling mirror 
in front of the laser, one can reflect this wasted light back in, allowing power 
to build up in the arms until the laser merely resupplies the mirror losses. 
This can dramatically reduce the power requirement for the laser. The first 
interferometers will work with laser powers of 5-10 W. This is attainable with 
modern laser technology. 

Quantum effects: Shot noise is a quantum noise, but in addition there 
are effects like bar detectors face: zero-point vibrations of mirror surfaces and 
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so on. These are small compared to present operating limits of detectors, but 
they may become important in 5 years or so. Practical schemes to reduce this 
noise have already been demonstrated in principle, but they need to be im- 
proved considerably. They can be reduced by making the mirror masses large, 
since the amplitude of vibration scales inversely as the square-root of the mass. 

Gravity gradient noise: One noise which cannot be screened out is that 
arises due to changes in the local Newtonian gravitational field on the time 
scale of the measurements. A gravitational wave detector will respond to tidal 
forces from local sources just as well as to gravitational waves. Environmental 
noise comes not only from man-made sources, but even more importantly from 
natural ones: seismic waves are accompanied by changes in the gravitational 
field, and changes in air pressure are accompanied by changes in air density. 
The spectrum falls steeply with increasing frequency, so for first-generation 
interferometers this will not be a problem, but it may limit the performance 
of detectors a decade from now. And it is the primary reason that detecting 
gravitational waves in the low-frequency band around 1 mHz must be done in 
space. 

2.4 Matched filtering and optimal signal-to- 
noise ratio 

Matched filtering is a data analysis technique that efficiently searches for a 
signal of known shape buried in noisy data (Helstrom, 1968). The technique 
consists in correlating the noisy output of each interferometer data with a set of 
theoretical waveform templates. If the maximum of the correlation exceeds a 
preset threshold, a detection is announced. However, even if a signal is present 
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in the data, the accompanying noise can mask its presence by preventing such 
crossing of the threshold in some cases. Thus a signal can be detected with 
only a certain probability, called its detection probability. Also, there will be 
a non-zero probability, called the false alarm probability, of a false detection 
due to noise alone. A formalism to choose the optimal set of templates us- 
ing a criterion was first given by Sathyaprakash and Dhurandhar (1991) and 
Dhurandhar and Sathyaprakash (1994) known as the S-D formalism. 

These search templates will be discrete subset of a continuous, multi pa- 
rameter family, each of which approximate a possible signal. Given a signal 
h(t) buried in noise n(t), the task is to find an 'optimal' template q(t) that 
would produce, on the average, the best possible SNR. The enhancement in 
the signal-to-noise ratio not only increases the number of detectable events 
but, more importantly, it also allows a more accurate determination of signal 
parameters — the error in the estimation of parameter being inversely propor- 
tional to the signal-to-noise ratio. 

If the functional form of the template is identical to that of the signal, the 
mean signal-to-noise ratio in the presence of a signal is highest possible for 
any linear data processing technique, which is why matched filtering is also 
known as optimal filtering (Helstrom, 1968). A matched filter , in the frequency 
domain, is a best-guess template of the expected signal waveform divided by 
the interferometer's spectral noise density. 

In order to take full advantage of matched filtering it is essential that 
the phase of the waveform be known to a very high degree of accuracy to 
accumulate the SNR (Cutler et al. 1993). This is referred to as coherent 
integration of the signal, where the SNR grows \/ 7 T . This is in contrast with 
incoherent methods, such as stacking and tracking (Brady, et. al. 1998, 2000), 
where one disregards the phase information and the SNR does not accumulate 
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as quickly. Tracking involves tracking of lines in the time-frequency plane 
built from the FT of short (around 40 minutes long) stretches of data. The 
other stacking involves dividing the data into shorter (around a day long) 
stretches, searching each stretch for signals, and enhancing the detectability 
by incoherently summing the FT of data stretches. 

There are several questions that must be answered in order to determine the 
feasibility of a matched filtering search strategy and, if feasible, to implement 
it. Which parameters significantly affect the wave form ? How should the 
spacing of the template parameters (lattice points) be chosen? Is there a 
parameterization that is in some sense preferred by the template waveforms? 
How many templates are needed to cover a given region of interest in the 
parameter space and how much computing power and memory will it cost to 
process the data through them? In the case of modulated Pulsar signal we 
must also ask what approximation to the true wave form is good enough? 

Ambiguity function, well known in statistical theory of signal detection 
(Helstrom, 1968), is a very powerful tool in signal analysis. It helps to access 
the number of templates required to span the parameter space of the signal. 
However, Apostolatos (1995) defined a term called Fitting Factor (FF) which 
quantitatively describes the closeness of the true signals to the template man- 
ifold in terms of the reduction in SNR due to cross correlation of a signal 
outside the manifold with all the templates lying inside the manifold. 

Fitting factor 

One of the standard measure for deciding what class of wave form is good 
enough is the Fitting Factor and has already been used by various earlier au- 
thors (Sathyaprakash and Dhurandhar, 1991 and Dhurandhar and Sathyaprakash, 



45 



1994) as a measure of adequateness of a template family, but previous discus- 
sions have not shown explicitly how this intuitively well-formed function is 
connected with detectability. Previously FF has sometimes been called the 
"correlation" , a number that is often used in so many different ways. The name 
comes from the fact that, by maximising the quantity FF, the parameters of 
the templates can be adjusted to best fit the true waveform. 

Let us consider a detector receiving almost monochromatic GW signals 
from a pulsar. The signal will be Doppler modulated due to the motion of 
Earth and stars, represented by S(t). The detector's output X(t) of the de- 
tector will be 

X(t) = S(t) + n(t) (2.7) 

If the search template is the exact waveform h(t), then we would achieve the 
highest possible signal-to-noise ratio given as 

S' 



v ; (S(f)\S(f)y (21 

where the inner product of two waveforms {hi\h 2 ) is defined as 

K(f)h2(f)+Hf)hl(f) 



(h\h 2 ) = 2 / 

JO 



S n (f) 



-df 



- «f» 

The SNR ratio will be reduced below (S/N) max whenever the template is 
not the exact waveform but some other approximate one. If search is performed 
with some family of templates S T (f;£j for a signal vector £ (The signal is 
considered a vector in parameter space). Then the SNR will be given by 

s\ = m f (x(f)\s T (f-Q) 

N) rms(n\S T (f;i)) 
For an ensemble of realizations of the detector noise, the expectation val- 
ues (n\S T (f;0) and (n\S T (f;g)){n\S T (f;g)) are zero and (S T (f;0\S T (f;0), 
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respectively. Thus the ensemble-averaged signal-to-noise ratio turns out to be 

■ s . ^ max {S (f)\S T (f-i)) 



Nl * ^(ST(f;i)\S T (f;i)) 

7 (S(f)\S T (f;£)) 



(2.11) 



(S T (f;0\S T (f;0)(S(f)\S(f))_ 
The name Fitting Factor is given to the reduction in signal-to-noise ratio that 
results because of using the chosen template family, rather than the true signal 
h, in the search: 

FF @ = T (S(f)\S T (f;0) 



^(S T (f;0\S T (f;0) 



= J (S(f)\S T (f;j)) 



(S T (f;0\SAf;0}(S(f)\S(f)} 
Therefore, if the family of search template used at the detection stage does 
not contain the true signal waveform, then the reduction in the signal-to-noise 
ratio will be given by 

f)^tt)x(|) (2.13) 

/ \ / max 

The spacing of the grid of filters is decided by the fractional loss due to to 
the imperfect match that can be tolerated. Given the parameters space that 
one needs to scan, it is then easy to estimate the total number of filters required 
to carry out the search for the signal. The signal strength scales inversely with 
distance, the fraction of event rate retained is approximately FF 3 . Therefore 
it has become conventional to regard FF = 0.9 as the lowest acceptable FF 
value for some template family to be considered adequate. Since a reduction 
in SNR by 10% means a 27% loss in the event rate. On the other hand a 10% 
reduction in SNR is equivalent to roughly 10% shortening of the detectors 
arms. 
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2.5 Computational costs 



Matched filtering places stringent demands on the knowledge of the signal's 
phase evolution which depends on two things: (1) our modelling of the signal 
and (2) the parameters characterising the signal. If our signal model is inaccu- 
rate or if the signal's parameters are unknown, there could be a loss in the SNR 
extracted. For instance, in the case of inspiral signals, a mismatch of one cycle 
in 10 4 cycles leads to a drop in the SNR by more than a factor of two, losing a 
factor of eight in the number of potentially detectable events. Recall that the 
SNR is inversely proportional to the distance to a source; thus an SNR loss 
by a factor a will reduce the span of a detector by the same factor, resulting 
in a decrease in the volume of observation, and hence the number of events, 
by a factor a 3 . Moreover, since the parameters of a signal will not be known 
in advance, it is necessary to filter the data with a family of templates located 
at various points in the parameter space e.g., placed on a lattice such that 
any signal will lie close enough to at least one of the templates to have a good 
cross-correlation with that template. The number of such templates is typi- 
cally very large. This places a great demand on the computational resources 
needed to make an on-line search. 

The search problem for continuous waves from spinning neutron stars is 
the most computer-intensive job in gravitational wave data analysis. Today, 
there is little hope that all-sky searches lasting for a year or more, can be 
made. It is easy to see why this is such an intensive job: Firstly, the data 
has to be collected continuously for months together and at a good sensitivity. 
Secondly, though a neutron star emits a periodic signal in its rest frame, save 
for the neutron star spin-down which indeed induces some modulation in the 
waveform, because of Earth's acceleration relative to the source, the detector 
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does not see a periodic wave. The wave is both frequency and amplitude mod- 
ulated. One can, fortunately, de-modulate these effects since Earth's motion 
is known quite accurately, and hence recover the original periodic signal. But 
de-modulation requires a knowledge of the source's direction and its frequency, 
which are unknown in a blind search. The angular resolution one obtains in a 
year's integration is A6 = X/D, where A is the wave length of radiation and D 
is the baseline of the detector in a year's integration, namely 1 A.U. Thus, for 
/ = 1 kHz we have A6 = 10~ 6 rad or about two arcsec. Now, assuming that 
the source may be in any one of the 4 arcsec 2 patches on the sky we get the 
number of patches in the sky for which we will have to try out a de-modulation 
correction to be 4n/(A8) 2 = 47rl0 12 . It is quite an impossible task to apply 
Doppler de-modulation to the detector output for each of these ~ 10 13 patches 
and compute as many Fourier transforms. 

One, therefore, asks the question given a compute power what is the best 
possible search one can do? Is there any advantage in going from a one-step 
search to a two or multi-step hierarchical search? What about directional 
searches? These are some of the problems for which we have some answer; but 
a great deal of work is needed and is currently under progress, to improve and 
optimise search algorithms. 

The differential geometric studies shows that with a TFLOPS-class com- 
puter, the number of days of data can be analysed on-line and we can carry 
out a blind search. Unfortunately, the longest data we can integrate on-line, 
for neutron stars with spin frequencies / < 100 Hz and spin-down rates less 
than 1000 years, is about 18 days. This yields a SNR lower by a factor of 
5 as compared to a year's worth of observing. On-line searches for neutron 
stars with / < 500 Hz (largest observed frequencies of millisecond pulsars) 
and spin-down rates of 40 years (shortest observed spin-down rates), can only 
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be made for a data set lasting for a duration of 20 hours or less. If source's 
position is known in advance, but not its frequency, then one can carry out an 
on-line search, again with a TFLOPS-class computer, for the frequency of the 
source in a data set that is worth 3 months long. This is good news since there 
are many known pulsars and X-ray binary systems that are potential sources 
of radiation. In addition, the obvious targeted search locations are the centre 
of the Galaxy and globular clusters. 

There have been efforts to study the effectualness of a two-step hierarchi- 
cal method for a blind search. Here the basic idea is to construct Fourier 
transforms of data sets of duration smaller than the period in which Doppler 
modulations will be important and to stack spectral densities obtained in this 
way and to add them all up. This is an incoherent way of building the signal 
since one adds spectral densities that have no phase information. Therefore, 
one gains in SNR less than what an optimal matched filtering method is able 
to achieve. However, this does not matter since (i) the targeted SNR's are 
quite high ~ 10 and (ii) candidate events can always be followed-up using 
coherent integration methods. These methods afford an on-line all-sky blind 
search for continuous gravitational waves for a period of 4 months or less using 
a 20 GFLOPS computer. Detector groups are planning to build computers of 
this kind to aid in their search for continuous GW. 

2.6 Detection criteria 

A signal needs to be above the noise experienced in the instrument and envi- 
ronment. However, this alone is insufficient to establish it as a gravitational 
wave in the terrestrial detectors. The most satisfying circumstance is that a 
gravitational wave observation be made in a set of widely distributed detectors 
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[the Gravitational Wave Network (GWN)] and the resorted wave forms allow 
the solutions for the polarization of the wave and the position of the source. 
Armed with this information an electromagnetic (or neutrino) search could be 
attempted in the error circle of the gravitational wave detection: a time hon- 
ored approach bringing gravitational wave of observation into the main stream 
of Astrophysics. The strategy would apply to all classes of sources: impulsive 
chirps, quasiperiodic and periodic. 

The confident detection of impulsive sources is most difficult, while the pe- 
riodic and quasiperiodic detection will have confidence limits based on quasi- 
stationary system noise (the signals last long enough to take a meaningful 
sample of the noise spectrum), the impulsive signals especially if rare, will 
be particularly dependent of the non-Gaussian component of the noise; the 
noise most difficult to reduce and control in a signal detectors. The technique 
of multiple coincidence of several detectors is one of the best means to gain 
confidence. The coincidence must occur within a time window to permit a 
consistent solution for a location in the sky. If the general character of the 
source can be guessed in advance (for example, a binary coalescence chirp, or 
a black- hole normal mode oscillation), the signal is filtered prior to the coinci- 
dence measurement to improve the sensitivity. The more detectors involved, 
the greater the confidence assigned to the detection. 

There is still a possibility of coincidence due to environmental or anthro- 
pogenic causes. The various sites throughout the world are far enough apart 
that most environmental perturbations should not correlate between them. 
The acoustic noise, the seismic noise, and the power line (especially if the 
network includes detectors in a different power grids and significantly differ- 
ent time zones) will be uncorrelated. There are correlations in the magnetic 
field fluctuations (thunderstorms) and in radio frequency emittions. As part 
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of the detection strategy a large number of environmental parameters will be 
measured along with the gravitational wave signals at each site. One of the 
requirements for the authenticity of impulsive sources will be the lack of cor- 
relation with environmental perturbations and other ancillary internal signals 
developed to monitor the performance of the instruments. 
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Chapter 3 



Data analysis of continuous 
gravitational wave: Fourier 
transform-I 

3.1 Introduction 

The first generation of long-baseline laser interferometers and ultra cryogenic 
bar detectors will start collecting data very soon. The network of detectors 
will not only confirm the existence of GWs but will also yield information 
about its important parameters viz; amplitude, frequency and phase. In ad- 
dition, a careful monitoring of the signal will also provide the information 
about the structure and the dynamics of its source. At the present stage, the 
data analysis depends largely upon the study of the expected characteristic 
of its potential sources and the waveforms. The majority of the experimental 
searches are focussed on the detection of burst and chirp signals. However, the 
interest in the data analysis for CGW signals is growing. A prime example of 
sources of this type is a spinning neutron star. Many research groups around 
the globe are working extensively on the data analysis for spinning neutron 
stars [Jaranwoski, et.al (1998, 1999, 2000), Brady et. al. (1998, 2000), Krolak 
(1999)]. 

Our Galaxy is expected to have, on logistic grounds, at least 10 8 spinning 
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neutron stars; with a birth rate of one in every 30 years. However, a very 
insignificant fraction of them are observed — the number of known pulsars has 
gone up to 1100. There are compelling arguments that nearby millisecond 
pulsars can provide a detectable source of CGW. Hence, there is a strong case 
for searching all sky for the presence of neutron stars with the philosophy that 
they are emitting GWs but otherwise unknown optically. This study will help 
us to understand the pulsar population in our Galaxy. 

The detection of GW signals in the output has its own problems, not 
the least of which is the sheer volume of data analysis. Bar detectors have 
essentially the same problems as interferometers in reference to CGW sources. 
The possibility that these detectors are capable to search for different kinds 
of signals makes them very versatile, but also complicates the data analysis. 
Each detector produces a single data stream that may contain many kinds 
of signals. Detectors don't point, but rather sweep their broad quadrupolar 
beam pattern across the sky as the earth moves. So possible sources could 
be anywhere on the sky and accordingly the data analysis algorithms need to 
accommodate signals from any arbitrary location of its source. 

In this and the next Chapter we present analysis of FT of the output data 
of a ground based laser interferometer. The output data has prominently 
dominant broad band noise and the signal is to be extracted out of it. For 
this, one has to enhance SNR. This is achieved by analyzing long observation 
time data as SNR is directly proportional to the square root of observation 
time \[T~ . However, in a data for long duration, the monochromatic signal 
gets Doppler modulated due to (i) orbital motions of Earth around Sun and 
(ii) spin of Earth. The Frequency modulation (FM) will spread the signal in a 
very large number of bins depending on the source location and the frequency. 
In addition there is Amplitude modulation (AM). As we will see in the sequel 
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the amplitude of the detector output consists of simple harmonic terms with 
frequencies w rot and 2w rot where, w rot stands for angular rotational frequency 
of Earth. Accordingly, the AM results in splitting of FT into frequencies ±w rot 
and ±2w rot . 

In the next section, we present the noise free response of the laser inter- 
ferometric detector and obtain the explicit beam pattern functions. In section 
3, we discuss the Doppler effect and obtain the FT of the FM signal for ar- 
bitrary source and detector locations taking into account Earth's rotational 
motion about its axis and its revolution around Sun. In section 4, the FT of 
the Doppler modulated complete response of the detector has been obtained. 
In the last section we discuss the results of the Chapter. 

3.2 The noise free response of detector: Beam 
pattern and amplitude modulation 

Let a plane GW falls on a laser interferometer and produces changes in the 
arms of the detector. In order to express these changes quantitatively we would 
require to specify the wave and the detector. Let XYZ and xyz represent 
respective frames characterising wave and detector. We assume the direction 
of propagation of the wave to be the Z axis and the vertical at place of detector 
to be the z axis [see Fig. (|3.1j) ]. In view of our discussion in subsection (jl.2.2|) . 
the difference of the changes 81 in the arm lengths of the detector is given via 

A/ 

R(t) = - = - sin 2tt \(A x x A y x - A Y A y Y ) h+ + (A X X A Y + A* A^) h x ) (3.1) 

'o 

where l D is the normal length of the arms of the detector and 2Q express the 
angle between them (Schutz and Tinto, 1987). The matrix (A^J represents 
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Figure 3.1: Orientation of the detector. 



Figure 3.2: Earth frame (x',y',z f ). 



Figure 3.3: Solar System Barycentre frame (X',Y',Z'). 
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the transformation expressing the rotations to bring the wave frame (X, Y, Z) 
to the detector frame (x, y, z). The direction of the source may be expressed 
in any of the coordinates employed in Spherical Astronomy. However, we find 
it convenient to define it in Solar System Barycentre (SSB) frame (X', Y', Z') 
[see Fig. (13. 3j) ]. This SSB frame is nothing but astronomer's ecliptic coordinate 
system. Let 9 and <f> denote the celestial colatitude and celestial longitude of 
the source. These coordinates are related to right ascension a and declination 
5 of the source via 

cos 9 = sin 5 cos e — cos S sin e sin a 
sin 9 cos (j) = cos 8 cos a > (3.2) 

sin 9 sin <\> = sin 5 sin e + cos 5 cos e sin a ) 



where e represents obliquity of the ecliptic. We choose x axis as the bisector 
of the angle between the arms of the detector. At this stage the orientation of 
the detector in the horizontal plane is arbitrary. It is assigned with the help of 
the angle 7 which x axis makes with the local meridian. The location of the 
detector on earth is characterised by the angle a, colatitude and (3, the local 
sidereal time expressed in radians. The transformation matrix (A^-J may be 
expressed as 

A = DCB (3.3) 

where 

B : rotation required to bring XYZ to X'Y'Z' 
C : rotation required to bring X'Y'Z' to x'y'z' 
D : rotation required to bring x'y'z' to xyz 

Here x'y'z' represents the frame associated with earth [see Figs. (j3.1H3.2|) ]. The 
Euler angles defining the corresponding rotation matrices (Goldstein, 1980) are 
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given via 



B : (0,0,V>) 1 

C : (0,e,0) 1 (3.4) 

D : (a,/3 + 7r/2, 7 -7r/2) J 

where ip is a measure of the polarisation of the wave. Let us express Eq. (j3.1|) 
as 

51 

R(t) = — = -sin2fi[F + /i+ + F x h x ] (3.5) 

The functions F + and F x involve the angles 9, 0, ip, e, a, j3, 7 and express 
the effect of the interaction of the wave and detector. These are called antenna 
or beam patterns. After a straight forward calculation one obtains (Jotania 
and Dhurandhar, 1994) 



F + (t) = - [{(cos-^cos^ — cos^sin^sin-?/') 2 — (cos-^sin^ + cos^cos^sin-?/') 2 } 
x {2 (cos a cos (3 cos 7 — sin (3 shi7)(— cos a cos (3 sin 7 — sin (3 cos 7)} 
+ {(— sin-0 cos — cos 9 sin (f)cosip) 2 — (cos# cos cos?/' — sin^ sin0) 2 j 
x [cos 2 e {2 (cos a sin (3 cos 7 + cos (3 sin 7) (cos (3 cos 7 — cos a sin (3 sin 7) } 

— sin 2 e sin 2 a sin 27 + sin 2e {(cos a sin (3 cos 7 + cos /3 sin 7) sin a sin 7 

— (cos (3 cos 7 — cos a sin /3 sin 7) sin a cos 7}] 

— (sin 2 9 cos 20) 

x [sin 2 e {2 (cos a sin f3 cos 7 + cos /3 sin 7) (cos f3 cos 7 — cos a sin (3 sin 7)} 

— cos 2 e sin 2 a sin 27 — sin 2e { (cos a sin (3 cos 7 + cos (3 sin 7) sin a sin 7 

— (cos {3 cos 7 — cos a sin /3 sin 7) sin a cos 7}]] 

+ {(cos^ cos0 — cos# sin 0sin-0)(— sin -0 cos — cos 9 sin cosip) 

— (cos ip sin + cos 9 cos sin ■0) (cos 9 cos cos ip — sin -0 sin 0) } 
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x [cos e { (cos (3 cos 7 — cos a sin (3 sin 7) (cos a cos (3 cos 7 — sin (3 sin 7) 

— (cos a sin /3 cos 7 + cos /3 sin 7) (cos a cos (3 sin 7 + sin /3 cos 7) } 
+ sin e { (cos a cos (3 cos 7 — sin (3 sin 7) (sin a sin 7) 

+ (cos a cos (3 sin 7 + sin /3 cos 7) (sin a cos 7) }] 

+ { (cos ip cos — cos 9 sin sin ijj) (sin # sin 0) 

+ (cos -0 sin + cos 6* cos sin ip) (sin 6* cos 0) } 

x [— sin e { (cos (3 cos 7 — cos a sin /3 sin 7) (cos a cos /3 cos 7 — sin /3 sin 7) 

— (cos a sin f3 cos 7 + cos (3 sin 7) (cos a cos (3 sin 7 + sin /3 cos 7) } 
+ cos e { (cos a cos /3 cos 7 — sin (3 sin 7) (sin a sin 7) 

+ (cos a cos /3 sin 7 + sin (3 cos 7) (sin a cos 7}] 

+ {(— sin -0 cos — cos 9 sin cos -0) (sin sin 0) 

+ (cos cos cos -0 — sin -0 sin 0) (sin 9 cos 0) } 

x [—(cose sine) {2 (cos a sin /5 cos 7 + cos (3 sin 7) 
(cos /3 cos 7 — cos a sin /3 sin 7) + sin 2 a sin 27 } 

+ cos 2e { (cos a sin (3 cos 7 + cos /3 sin 7) (sin a sin 7) 

— (cos (3 cos 7 — cos a sin /3 sin 7) (sin a cos 7)}] ; (3.6) 

F x (t) = 2(cos-0 cos0 — cos 0sin sin-0)(cos-0 sin + cos0 cos0 sin-0) 

x (cos a cos p cos 7 — sin (3 sin 7) (— cos a cos (3 sin 7 — sin /3 cos 7) 
+ (— sin -0 cos — cos 9 sin cos -0) (cos 6* cos cos -0 — sin ip sin 0) 
x [cos 2 e {2 (cos a sin (3 cos 7 + cos (3 sin 7) (cos /3 cos 7 — cos a sin f3 sin 7)} 

— sin 2 e sin 2 a sin 27 + sin 2e { (cos a sin (3 cos 7 + cos /3 sin 7) sin a sin 7 

— (cos (3 cos 7 — cos a sin /5 sin 7) sin a cos 7}] 

— - (sin 2 9 sin 20) 

59 



x [sin 2 e {2 (cos a sin (3 cos 7 + cos (3 sin 7) (cos (3 cos 7 — cos a sin (3 sin 7)} 

— cos 2 e sin 2 a sin 27 — sin 2e { (cos a sin (3 cos 7 + cos /? sin 7) sin a sin 7 

— (cos (3 cos 7 — cos a sin /3 sin 7) sin a cos 7}] 

+ {(cos-0 cos0 — COS6 1 sin 0sin-0)(— sin ^ sin + cos 9 cos cos -0) 

— (cos -0 sin + cos 9 cos sin ip) (cos 6* sin cos -0 + sin -0 cos 0) } 

x [cos e { (cos /3 cos 7 — cos a sin (3 sin 7) (cos a cos (3 cos 7 — sin f3 sin 7) 

— (cos a sin /3 cos 7 + cos /3 sin 7) (cos a cos /3 sin 7 + sin /? cos 7) } 
+ sin e { (cos a cos /3 cos 7 — sin /3 sin 7) (sin a sin 7) 

+ (cos a cos (3 sin 7 + sin /3 cos 7) (sin a cos 7) }] 
+ { — (cos -0 cos — cos 9 sin sin -0) (sin # cos 0) 
+ (cos -0 cos + cos 9 cos sin -0) (sin 6* sin 0) } 

x [— sin e { (cos (3 cos 7 — cos a sin (3 sin 7) (cos a cos (3 cos 7 — sin f3 sin 7) 

— (cos a sin /3 cos 7 + cos /3 sin 7) (cos a cos /3 sin 7 + sin f3 cos 7) } 
+ cos e { (cos a cos /3 cos 7 — sin (3 sin 7) (sin a sin 7) 

+ (cos a cos /3 sin 7 + sin /3 cos 7) (sin a cos 7}] 
+ {(— sin -0 sin + cos 6* cos 0cos-0)(sin# sin 0) 
+ (cos 9 sin cos -0 + sin -0 cos 0+) (sin 9 cos 0)} 
x [— (cos e sin e) {2 (cos a sin /3 cos 7 + cos /3 sin 7) 
(cos (3 cos 7 — cos a: sin (3 sin 7) + sin 2 a sin 27 } 
+ cos 2e { (cos a sin f3 cos 7 + cos (3 sin 7) (sin a sin 7) 

— (cos /9 cos 7 — cos a sin /3 sin 7) (sin a: cos 7)}] ; (3.7) 

It is easy to understand the structure of the above complicated looking ex- 
pressions for F + and F x by introducing following abbreviations 
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U = cos a cos (3 cos 7 — sin (5 sin 7 , 

V = — cos a cos (3 sin 7 — sin (5 cos 7 , 

X = cos a sin /3 cos 7 + cos /3 sin 7 , 

F = — cos a sin /3 sin 7 + cos /3 cos 7 



(3.8) 



L = cos ip cos — cos 9 sin sin tp , 

M = cos ip sin + cos^ cos sin ^ , 

N = — sin -0 cos — COS6 1 sin cos -0 , 

P = — sin -0 sin + cos9 cos cos -0 , 

Q = sin 9 sin , R — sin cos , 



(3.9) 



A = 

B = 

C = 

D = 

E = 



2XY cos 2 e — sin 2 e sin 2 a sin 27 + sin 2e(X sin a sin 7 

— F sin a cos 7) , 

2XF sin 2 e — cos 2 e sin 2 a sin 27 — sin 2e(X sin a sin 7 

— F sin a cos 7) , 

cos e(YU + XV) + sin e(U sin a sin 7 — V sin a cos 7) , 

— sin e(YU + XV) + cos e(£7 sin a sin 7 — V sin a cos 7] , 
— 2XY cos e sin e — cos e sin e sin 2 a sin 27 

+ cos 2e(X sin a sin 7 — F sin a cos 7) 



> (3.10) 



One obtains: 



F + (t) 



- [2{L 2 - M 2 )f/F + {N 2 - P 2 )A + (Q 2 - P 2 )P] + (LJV - MP)C 



+ (LQ + MR)D + (JVQ + PR)E , 



F x (t) = 2LMUV + iVPA - -P sin 2 9 sin 20 + (LP + MiV)C 
+ (MQ - LP)P + (Pg - NR)E 



(3.11) 



(3.12) 



The compactification achieved here is because of the fact that these abbrevi- 
ations find places in the transformation matrices as: 
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D 



L N Q 

B | M P -R 

sin 6 sin ip sin 9 cos ip cos # 

1 

C = I cose sine 

— sin e cos e 

U V sin a cos (3 

X Y sin a sin (3 

sin a cos 7 sin a sin 7 cos a 



(3.13) 



(3.14) 



(3.15) 



After algebraic manipulation Eqs. ()3.11|) and ()3.12j) may be expressed as 

F+(t) = F 1+ cos 2(3 + F 2+ sin 2/3 + F 3+ cos /3 + F 4+ sin /3 + F 5+ ; (3.16) 

F x (t) = F lx cos2(3 + F 2x sin2/? + F 3x cos/3 + F 4x sin/3 + F 5x (3.17) 
where F, + and F ix (i = 1, 2, 3, 4, 5) are time independent expressions given via 



2+ 



'3+ 



-2G cos a cos 2 7 + ^f^(cos 2 a + 1) , ^ 
H cos a cos 27 + G sin 27(cos 2 a + 1) , 
J sin a cos 27 + J sin 2a sin 27 , 
2 J sin a cos 27 — | sin 2a sin 27 , 



G = \ [{LQ + MR) sin e - (ZJV - MP) cos e] , 

# = I[(Ar2_p2 )cos 2 e _ (L 2_ M 2 ) + ( g2_ jR 2 )sin 2 e 

-(iVQ + PP) sin2e], 

/ = i[(Q2 - J R 2 ) s in2e- (iV 2 -P 2 )sin2e- 

2(NQ + PR) cos2e], 
J = \[{LN — MP) sin e + (LQ + MR) cos e] 



(3.18) 



(3.19) 



Let us note that F ix is related to F i+ via 



F ix {9,(f),ip,a,^,^,e) 



71 



(3.20) 



1,2,3,4,5 
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This symmetry is representative of the quadrupolar nature of the detector and 
the wave. A detector in different orientations will record different amplitudes 
in the response. The explicit beam pattern functions may be computed easily 
for any instant of time. Due to symmetries involved in F + and F x it is sufficient 
to discuss either of the beam pattern. 

The amplitude modulation of the received signal is a direct consequence 
of the non-uniformity of the sensitivity pattern. The antenna response is not 
peaked in any direction but is anisotropic. As remarked earlier they are fairly 
complicated function of their arguments. Equations. ()3.16|) and ()3.17|) reveal 
that the monochromatic signal frequency will split, due to AM, into five lines. 
This results in the distribution of energy in various frequencies and consequent 
reduction of the amplitude of the signal. The periodicity of the beam patterns 
F + and F x with a period equal to one sidereal day is due to the diurnal motion 
of Earth. 

3.3 Doppler shift and Frequency modulation 

The frequency of a monochromatic signal will be Doppler shifted due to the 
translatory motion of the detector, acquired from the motions of Earth. Let 
us consider a CG W signal of constant frequency f Q . The frequency /' received 
at the instant t by the detector is given by 



where n is the unit vector from the antenna to source, v is the relative velocity 
of the source and the antenna, and c is the velocity of light. The unit vector 
n from the antenna to the source, because of the fact that the distance of the 




(3.21) 
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source is very large compared to the average distance of the centre of the SSB 
frame and the detector, may be taken parallel to the unit vector drawn from 
the centre of the SSB frame to the source. Hence, 



n = (sin 9 cos , sin 9 sin , cos#) (3.22) 

As v keeps on changing continuously both in its amplitude and direction /' is 
a continuous function of t. Further, since t;<cwe take j — l. 
The radius vector r(t) in the SSB frame is given by 

r(t) = [R se cos(w or bt) + R e sin a cos/? , 

R se sm.{w or bt) + R e sin a sin (3 cos e — R e cos a sin e , 

R e sin a sin (3 sin e + R e cos a cos e] ; (3.23) 

/3 = O + w rot t (3.24) 

where R e , R se and w or 6 represent respectively Earth's radius, average distance 

between Earth's centre from the centre of SSB frame and the orbital angular 

velocity of Earth. Here, t represents the time in seconds elapsed from the 

instant Sun is at the Vernal Equinox and f3 a is local sidereal time at that 

instant. The Doppler shift is now given via 

f'-fo v.n r t .n R se w orb . . 

— -„ = it) = = sin0sin(0 - w orb t) + 

Jo c c c 

eWrot s in a [sin ^{cos (3 cos e sin 0— 

c 

cos sin (3} + cos (3 sine cos 9] (3.25) 

The phase of the received signal is given by 

$(f) = 2tt [ f'(t')dt' 
Jo 

= 2nf \l + —(If)] dt' (3.26) 



o 
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Here, we assume the initial phase of the wave to be zero. After a straight 
forward calculation we obtain 



27Tfo 



t + sin 9 cos <fi'+ 
c 



R P 



sin a{sin 9 (sin (3 cos e sin + cos cos (3) + 

Rf, 



sin (3 sin e cos 9} sin 9 cos — 

c 

R 

— sin a{sin #(sin (3 a cos e sin + cos cos (3 Q ) + 
sin /9 sine cos 9}} 

Zftfot + Z COs(w or bt — 0) + V Sin(w rot t) + Q COs(w rot t) — 1Z — Q 

2rcf t + Z cos« rot - 0) + M cos(£ rot -6)-K-Q (3.27) 



where 



Q 

N 
z 

K 

5 

£,orb 



27r/ ^ sin a (cos /3 (sin cos e sin + cos 9 sin e) 
— sin f3 sin 9 cos 0) , 

27r/ -^ sin a (sin /3 G (sin cos e sin + cos 9 sin e) 
+ cos f3 Q sin cos 0) , 

V? 72 + Q 2 , 
27r/ o ^sin0, 

Z COS , 



tan 



-1 v 



w orb t - , 

W or bt = d^rot't a 
W ro tt 



w orb /w. 



rot 



1/365.26, 



(3.28) 



(3.29) 



The two polarisation states of the signal can be taken as 



h + (t) = h 0+ cos[$(t)] 



h x (t)=h 0x sin[$(*)]; 



(3.30) 
(3.31) 
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h 0+ , h 0x are the time independent amplitude of h + (t), and h x (t) respectively. 

To understand the nature of the FM, let us consider the FT of h + (t) and 
h x (t). We may begin by considering the function 



h(t) = cos [$(*)] 



(3.32) 



and analysing the data for one day observation time. The FT is given via 



kf) 



cos[$(t)}e- i27rft dt ■ 



(3.33) 



T = one sideral day = 86164 sec. 
This may be splitted into two terms as 



Kf) 



(3.34) 



2w 



rot JO 



2tt 



3 i[fi/_+.Zcos(a£-<£)+A/ r cos(£-<S)-7e.-g] 



2w. 

fo±f 



rot 



(3.35) 



— / ^ e -i[^++Zcos(a^)+Afcos^-S)-TZ~Q] d ^ (3.36) 
'Jrot Jo 

€ = irot = w rot t (3.37) 



Numerical result shows that I v+ oscillates very fast and contributes very little 



to 



h(f) . Hence, hereafter, we drop /„, from Eq. ()3.34j) and write v in place 



of Using the identity 



i=oo 



„±jft COS 1? 



J (±k) + 2^Vjj(±/c) cos/0 



(3.38) 



we obtain 
Kf) . 



2m 



2vr 

_ e i{-n-Q) I e ivt 



rot 



'0 

m=oo 



k=oo 



J {Z) + 2 J k{Z)i k cos fc(a£ 



k=i 



J {N) + 2 Jm{M)i m cosm(£ - 5) 



m=l 



(3.39) 
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where J stands for the Bessel function of first kind of integral order. After 
performing the integration we get 



A 
B 
C 

V 



Mf) 



v 



k=oo m=oo 



2m 



rot 



E E - iV \ ; 



k=— oo m=—oo 



(k+m)ir 



K-Q 



J k (Z).J m (N) 

v 2 — (afc+m) 2 

sin 2vk cos(2akir — kef) — m5) — 

afc+m j CQS 2 U7T s in(2aA;7r — kef) — mS) + sm(k(f> + m5) } 

cos 2i/7r cos(2afc7r — k<f> — mS) + 

ka+m sin 2i/7r sin(2aA;7r — k<f> — mS) — cos(kcf) + mS) 



(3.40) 



(3.41) 



The FT of the two polarisation states of the wave can now be written as 



M/) 



vh 



fc=oo m=oo 



°+ 



2w 



rot 



k=—oo m=—oo 



h(f)} = -ih Qx \h(f)} 

- d L J ( 



fc=oo m=oo 



2w r , 



"EE e M i3[P-^C] 

rot , 



fc=— oo m=— oo 



(3.42) 



(3.43) 



The FT of the FM signal contains the double series Bessel functions. The 
Bessel functions has contributions due to the rotational as well as the orbital 
motion of the earth. It is remarked that Jotania et al. (1996) have analysed FT 
of FM signal for one day observation time. They have taken specific detector 
as well as source location. They have also neglected the orbital motion. Our 
analysis generalizes their results. We may now compute h(f) and may plot 

J d 

its behaviour. To achieve this we have made use of Mathematica software. 
Fig. (j3.4J) represents such a plot for 
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f a = 8QHz, h a = h x = l } 

a = n/3, /3 = tt/4, 7 = 2tt/5, > (3.44) 

9 = n/36, 4> = 7r , ip — 7t/6. ) 

with the resolution equal to l/T Q = 1.16 x 10~ 5 Hz. Figures (|H.5jl and ()3.6j) 
represent the plot of the FT at resolution 10~ 6 Hz and 10~ 7 Hz. A careful 
look at these plots reveals that the resolution of Fig. (|H.4|) does not represent 
the details of dominant peaks around f a , whereas, Fig. (|H.6jl does not give any- 
new behaviour as compared to Fig. (|H.5jl . Hence, we may say that a resolution 
of about 10~ 6 Hz is required to understand the correct behaviour of the FT 
for one day observation time data analysis. In this reference let us recall 
that the data analysis for Fast Fourier Transform (FFT) limits the resolution 
to 1/T Q . However, detector output may provide us higher resolution. Thus 
the semi-analytical analysis presented here may provide more information as 
compared to FFT. Further, in order to see the distribution of amplitudes at 
various frequencies we plot power spectrum in Fig. (|3.7|) . 
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Figure 3.4: FT of a FM signal of a source located at (n/36, n) with a resolution 
of 1.16 x 1CT 5 . 
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Figure 3.5: FT of a FM signal of a source located at (vr/36, ir) with a resolution 
of 10~ 6 . 
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Figure 3.6: FT of a FM signal of a source located at (7r/36, 7r) with a resolution 
of 10" 7 . 
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3.4 Fourier transform of the complete response 

The complete response R(t), in view of Eqs. (jSHJ), (pHKJ) . flTT7|) . (|3~30l and (ICTJ) 

may be written as 

R(t) = R+(t) + R x (t) ; (3.45) 

7? + (t) = /i Q+ [F 1+ cos 2/5 + F 2+ sin 2(3 + F 3+ cos (3+ 

F 4+ sin/? + F 5+ ] cos[$(f)] , (3.46) 
(*) = h 0x [F lx cos 2(3 + F 2x sin 2/? + F 3x cos/?+ 

F 4x sin/5 + F 5x ]sin[$(t)] (3.47) 

Here we have taken for simplicity the angles between the arms of the detector 
to be tt/2 i.e. Q = it /A. Now the FT of the complete response may be 
expressed as 



R(f) = R + (f)+R x (f) 

To evaluate these let us substitute (3 as given by ()3.24j) . One obtains 
'1 



(3.48) 



R + (t) 



"0+ 



{e- i2f3 °(F l+ + iF 2+ )e- i2Wrott + e i2f3 °(F l+ - iF 2+ )e i2Wrott + 



e~ iPo (F 3+ + iF 4+ )e- iWrott + e i(3 °(F 3+ - tF 4+ )e iWrott } 
+ F 5+ ] cos[$(t)] 



(3.49) 



and similar expression for R x {t). Now it is straight forward to obtain the 
expression for R+(f) and R x (f). We get 



R + (f) 



h 



e-W°(F 1+ +iF 2+ ) h(f + 2f rot )/2 



°+ 



J d 



+ 



h(f - 2f TOt )/2 



e— (F 1+ -iF 2 
e -i0o( F3++iFi+) ^(/ + / TOt )/2 

e^(F 3+ -zF 4+ ) \h(f-f rot )/2 



J d 



+ 



+ 



+ F, 



Kf) (3-50) 

d- 
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Similarly we have 



R*(f) 



h„ 



-iWo i 



(F 2x -iF lx ) h(f + 2f rot )/2 



+ 



_ e Ofc(j? +iFlx ) h(f-2f rot )/2 



+ 



{Fa 



iF 3 



+ / roi )/2 



e ^(F 4x +zF 3x ) h(f-f rot )/2 



iF R 



h(f) (3.51) 



Collecting our results the FT of the complete response of the detector for one 
day integration will be 



R(f) 



J2/3 



-i/3 



h(f + 2f rot )/2 [h 0+ (F 1+ +iF 2+ ) + h 0x (F 2x -iF 1> 

. d 

h(f - 2f rot )/2] [h 0+ (F 1+ - iF 2+ ) - h Qx (F 2x + iF lx )] + 

- d 

h(f + f rot )/2\ [h 0+ (F 3+ +iF 4+ ) + h 0x (F 4x -tF 3x )} + 

. d 

Hf - frot)/2\ j [K + ( F H - - h 0x (F 4x + tF 3x )} + 

(3.52) 



J d 



Mf) \h 0+ F 5+ - ih 0x F 5x ] 

. a 



This shows that due to AM every Doppler modulated FM signal will split 
in four additional lines at / ± 2f rot and / ± f ro t, where f rot is the rotational 
frequency of Earth (f ro t ~ 116 x 10~ 5 Hz). In terms of dimensionless scaled 
frequency v the power of the signal will be at v ± 2a, v ± a Sz v. 

We have plotted in Fig. ()3.8j) the power spectrum of the noise free complete 
response of the signal for its various parameters as given by 1)3.44)1 . The contri- 
bution in the power spectrum of modulation at frequencies / + 2f rot , f — 2f rot , 
f + f ro t and / - f rot and / are represented in Figs. (JH3|, IjTTTty . jSHJ), ((Hill 
and ()3.13j) . It is observed that the most of the power will be at / + 2f rot and 
least power will be in / — f rot . 
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Figure 3.7: Power spectrum of a FM signal of a source located at (7r/36,7r) 
with a resolution of 1CT 7 . 
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Figure 3.8: Power spectrum of the complete response of a Doppler modulated 
signal of a source located at (7r/36,7r) with a resolution of 1CT 7 . 
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Figure 3.9: Power spectrum of a Doppler modulated signal at frequencies 
/ + If rot of a source located at (n/36, ir) with a resolution of 1CT 7 . 
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Figure 3.10: Power spectrum of a Doppler modulated signal at frequencies 
/ — 2 f rot of a source located at (n/36, ir) with a resolution of 1CT 7 . 
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Figure 3.11: Power spectrum of a Doppler modulated signal at frequencies 
/ + frot of a source located at (n/36, it) with a resolution of 1CT 7 . 
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gure 3.12: Power spectrum of a Doppler modulated signal at frequencies 
— frot of a source located at (n/36, ir) with a resolution of 1CT 7 . 
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Figure 3.13: Power spectrum of a Doppler modulated signal at frequencies / 
of a source located at (vr/36, ir) with a resolution of 10 -7 . 
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3.5 Discussion 

In this Chapter we have considered the effect of Earth's motion on the response 
of the detector through FT analysis. It can easily be inferred from Eqs. (|3.52|) 
and (|3.421 13.43|) that the splitting of frequencies (i) in AM arise explicitly due 
to rotational motion and (ii) in FM arise due to rotational as well as orbital 
motion of earth. In view of the fact that the data output at the detector is 
available in discrete form, the analytical FT is not very convenient and one 
normally employs the popular FFT. However, FFT has resolution limited to 
1/T . Further, it is important to understand for how much time one can ignore 
the frequency shift arising due to Doppler effect. In fact, Schutz (1991) has 
demonstrated that these effects due to rotational motion are important after 
the time given by 



This means that for GW signal for frequency 80 Hz one has to take into 
account these effects after data time ~ 4 hours. The analytical FT studied in 
this Chapter leads to following inferences: 

(i) FFT for one day observation data will not provide sufficient resolution 
as to represent the correct picture of the frequency splitting. 

(ii) The adequate resolution required for one day observation ~ 10~ 6 Hz. 

(hi) The frequency split due to FM for frequency f a = 80 Hz and source at 
(0, 0) = (vr/36, 7r) is ~ 2 x 10~ 4 Hz and due to AM is ~ 4.64 x 10~ 5 Hz. 

(iv) The drop in amplitude due to FM alone is about 56% 

(v) The drop in amplitude due to AM alone is about 18%. 




(3.53) 
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(vi) The drop in amplitude for the complete response is about 74%. 

(vii) The maximum power due to AM is associated with f + 2f rot . 

It is remarked that the drop of the amplitude in complete response is severe 
both due to AM and FM as the relevant frequency range lie in the same region 
[see conclusion (iii) above]. 

Finally, we would like to mention that we have presented the FT analysis 
assuming the phase of the GW to be zero at that instant t — 0. However, 
one may relax this condition and may obtain the results easily by taking into 
consideration the effects of change of the time origin. 
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Chapter 4 



Data analysis of continuous 
gravitational wave: Fourier 
transform-II 

4.1 Introduction 

In previous Chapter we have presented Fourier analysis of one day observation 
data set of the response of a laser interferometer. We have seen that the 
amplitude and frequency modulations result into a large number of side bands 
about the signal frequency f Q . Consequently the maximum power lies in the 
frequency / + 2f rot with amplitude reduction by 74% to what one would have 
expected due to increased data interval. Hence, for GW detection it is desirable 
to obtain FT for larger data observation time interval. To begin with we 
present in the next section an analysis for data set of one year observation 
time. The analysis turns out to be simple in view of the observation that 
there exists correspondences and identifications to the analysis of previous 
Chapter to make them identical. We obtain in section 3 generalisations of the 
results for arbitrary observation time. To facilitate analogous modifications 
we have introduced corresponding quantities with tilde viz., C and V in place 
of C and T>. As an application of the results obtained here we consider spin 
down and N-component signal analysis respectively in sections 4 and 5. We 
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present discussion and conclusions in the last section. 



4.2 Fourier transform for one year integration 
4.2.1 Frequency modulation 

Let us rewrite the expression for phase of the GW signal referring to Eqs. (|3.27J1 . 
(E| and (El 

$(t) = 2nf t + Z cos(w arb t - 0) +Afcos(w rot t - 5) - 11 - Q (4.1) 
The FT for one year observation time T b s is given as 



h(f) = / cosmt)}e- l2lTft dt ; 
iy Jo 

a = a -1 = w rot /w or b ; T = one sideral day 
T nhq = aT ~ 3.14 x 10 7 sec. 



This may be splitted as usual into two terms as 



Kf) 



l/- I *V + ) 



j _ 1 /" e i[5P_+2cos(5-0)+ATcos(a5-5)--^-Q]^ 
2u>or& Jo 



1 



2tt 



2u>orfc JO 

/o±/ 
forb 



-i[iv++Z cos(£-<t>)+N 'cos(a£-<S)-72.-Q] ^ 



£ — £or& — Worbt 



(4.2) 
(4.3) 
(4.4) 



(4.5) 

(4.6) 
(4.7) 
(4.8) 



Hereafter, we neglect the contribution of I p+ to h(f) as it oscillates rapidly 

L J y 

and contributes very little and write v in place of P_. A careful comparison of 
Eq. (|4.6|) with (|3.35|) reveals that the integrand of the equations are identical 
with following identifications and correspondences. 
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5 ^ (j) , 
z <-> N 



(4.9) 



a. 



Hence, we may employ the results obtained there by introducing obvious cor- 
responding quantities i.e. B, C, V in place of B, C, V leaving A unchanged. 
Thus 



k=oo m=oo 



k=—oo m=— oo 



(4.10) 



A = 

B = 

C = 

V = 



J k (N)Jm(Z) 

u 2 -(ak+m) 2 ' 

sin2z/7rcos(2aA;7r — k5 — mtp) — 

^±^{cos 2z/7r sin(2a/c7r — k5 — m<j>) + sin(/c5 + m<j))} , 
cos2z/7rcos(2a/c7r — kS — m(f>)+ 

sin 2uix sin(2aA;7r — k5 — mcf)) — cos(/c5 + m0) 



) (4.H) 



ka+m 
v 



Now the FT of the two polarisation states can be written as 



and 



K(f) 



h 



J v 



fc=oo m=oo 



2w 



orb 



E E 



fc=— oo m=— oo 



-ifcox [*(/)] y 



2w orb 



k=oo m=oo 



E E *"®p-A 



k=—oo m=—oo 



(4.12) 



(4.13) 



The FT obtained contains double series of Bessel functions of the order k 
and m ranging from — oo to oo. It is well known that Bessel functions decrease 
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rapidly as the order exceeds the argument. Hence possible range of k and m 
over which the summation of the series is to be considered depends on the 
arguments of Bessel functions i.e Z and M . Referring to Eq. (|3.28J1 it is found 
that 



Zmax — 3133215 ( ^ KHz ) 

M max = 134 (dfc) 



(4.14) 



The FT of a FM signal for 



fo = 50 Hz , h = h x = 1 

a = 7r/4, P o = 0, 7 = vr, } (l.LV) 

= 7r/18, = 0, ^ = vr/4. 

is shown in Fig. (|4.1)1 . The spectrum has a resolution of 1/T D w 3.17 x 10^ 8 i?2;. 
The corresponding power spectrum is shown in Fig. (j4.2jl . We have convinced 
ourselves by plotting the FT at higher resolutions that resolution of 1/T G is 
sufficient to represent relevant peaks. We notice that the drop in amplitude is 
about 98%. This may be attributed to the presence of a very large number of 
side bands. 
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4.2.2 Complete response 



The complete response of the detector may be obtained employing Eqs. (j4.12j) . 
(ETT31) . (EraHJ) and (ITTTj) and one gets 

e -W°( Fl+ +iF 2+ ) \h(f + 2f rot )/2 



R + (f) 



h. 



°+ 



+ 



h(f - 2f rot )/2 



e iW °{F l+ -iF 2 
e -^(F 3+ +zF 4+ ) [h(f + f rot )/2 

e^°{F 3+ -iF A+ ) \h(f-f rot )/2 



+ 



+ 



+ F 5 



Hf) 



;(4.16) 



£x(/) 



e-^^-iFxJ h(f + 2f rot )/2 



J y 



e"^(F 4x -*F 3x ) [M/ + / rot )/2 
e ^(F 4x +iF 3 J [^(/-/ rot )/2 
After rearranging these expressions we obtain 



+ 



J y 



Hf) 



J y 



(4.17) 



*(/) 



-i2(3 



J y 



~h(f + 2f rot )/2 [ ho+ (F 1+ +iF 2+ ) + h 0x (F 2x -iF lx )] + 
J y 

h{f - 2f rot )/2\ [h 0+ (F 1+ - iF 2+ ) - h 0y (F 2x + iF lx )] + 
J y 

Hf + frot)/2] [h 0+ (F 3+ +iF A+ )+h 0x (F 4x -iF 3x )] + 
J y 

Hf - frot)/2\ [K + (F 3+ - iF i+ ) - h 0x {F ix + iF 3x )] + 
J y 

Hf)} [h 0+ F, + -ih 0x F, x ] (4.18) 

Figure ()4.3|) shows the power spectrum of the complete response of the 
Doppler modulated signal. We have kept here all parameters same as in FM. 
Figures ()4.4|) . ()4.5|) . ()4.6|) . ()4.7j) . and f)4.8|) are the power spectrum of complete 
response at frequencies / + 2f rot , f - 2f rot , f + f rot , f - f rot and /. In this 
case also we observe that the maximum power is associated with / + 2f rot and 
the least lies in / — f rot . 
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Figure 4.4: Power spectrum of a Doppler modulated signal at frequencies 
/ + 2f rot of a source located at (n/18, 0) with a resolution of 3.17 x 1CT 8 . 
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Figure 4.5: Power spectrum of a Doppler modulated signal at frequencies 
/ — 2f rot of a source located at (n/lS, 0) with a resolution of 3.17 x 10~ 8 . 
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Figure 4.7: Power spectrum of a Doppler modulated signal at frequencies 
/ — frot of a source located at (n/18, 0) with a resolution of 3.17 x 10~ 8 . 
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Figure 4.8: Power spectrum of a Doppler modulated signal at frequencies / of 
a source located at (7r/18, 0) with a resolution of 3.17 x 10~ 8 . 
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4.3 Fourier transform for arbitrary observa- 
tion time 



It is important to obtain the FT for arbitrary observation time. In this section 
we present this analysis. The results obtained will be employed to outline how 
spin down of a pulsar due to gravitational radiation back reaction or due to 
some other mechanism can be taken into account. 
The FT for data of observation time T Q is given via 



h(f) = [ ° cos[$(t)]e~ i27Tft dt 
Jo 



(4.19) 



As usual, this splits into two terms as 



h{f) = I v _ + I v+ ■ 



i[&-+Zoos(a£-<l>)+Mcos(ji-5)-n-Q\ 



2w ro t Jq 

1 



-i[£v + +Z cos(a£-0)+.A/'cos(£— S)—Tl- 



2lV ro t Jo 
fo±f . 
frot 



(4.20) 

(4.21) 
(4.22) 
(4.23) 



As I v+ contributes very little to h(f) we drop I v+ and write u in place of z/_. 
Using the identity ()3.38|) we get 

fc=oo 



Kf) 



1 



'0 

m=oo 



k=l 



2w rot 

J {M) + 2 J m W m cosm{Z - S) 



J {Z) + 2 Jk(Z)i k cos k(a$, 



m=l 



(4.24) 



After performing the integration and proceeding in a straight forward manner 
we have 
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fc=oo m=oo 

h w - 9^- E E e ^ c ~ - ^ ; ( 4 - 25 ) 

k=—oo m=—oo 



2 

J k (Z)J m (N) 
v 2 — (ak+m) 2 



B 

C = sin z/£ cos(a/c£ + m£ G — k<p — m5) — 

afc+m | CQS ^ sin (afc^o + m£ G — k<\> — m8) + sin(A;0 + mS)} 
V = cos u^ cos(a/c£ + m£ c — k<p — m5) + 

ka+m sin i/£ sin(a/c£ + m£ c — k<p — mS) — cos(A;0 + m5) 

(4.26) 

The FT of the two polarisation states of the wave can now be written as 
h+(f) = h 0+ h(f) 

, k=oo m=oo 

= ^EE e^B\&-iV\; (4.27) 

k=—oo m=~oo 

K{f) = -ihoMf) 

, k=oo m=oo 

E e ^ - ^ ( 4 - 28 ) 

k=—oo m=—oo 

Now it is simple matter to obtain the FT of complete response. One gets 

R(f) = e- i2 ?°h(f + 2f rot )/2 [h 0+ (F 1+ +iF 2+ )+h 0x (F 2x - iF lx )] + 
e i2 ?°h(f - 2f rot )/2 [h 0+ (F 1+ - iF 2+ ) - h 0x (F 2x + iF lx )] + 
e- l H(f + f rot )/2 [h 0+ (F 3+ +iF 4+ ) + h 0x (F 4x -iF 3x )] + 
e*°h(f - frot)/2 [h Q+ (F 3+ - iF A+ ) - h 0x (F 4x + iF 3x )] + 
Hf) [K + F 5+ -th 0x F 5x ] (4.29) 



97 



j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 



T o = 120 days 




f -> 



Figure 4.9: FT of a FM signal of a source located at (vr/9, 7r/4) with a resolu- 
tion of 9.67 x 10 -8 . 
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Figure 4.10: FT of the complete response of a Doppler modulated of a source 
located at (n/9, 7r/4) with a resolution of 9.67 x 10~ 8 . 
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Figure 4.12: Power spectrum of the complete response of a Doppler modulated 
signal of a source located at (n/9, 7r/4) with a resolution of 9.67 x 10" 8 . 
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Figure 4.13: Power spectrum of a Doppler modulated signal at frequencies 
/ + If rot of a source located at (n/9, 7r/4) with a resolution of 9.67 x 10~ 8 . 
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Figure 4.14: Power spectrum of a Doppler modulated signal at frequencies 
/ — 2 f rot of a source located at (n/9, 7r/4) with a resolution of 9.67 x 10~ 8 . 
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Figure 4.16: Power spectrum of a Doppler modulated signal at frequencies 
/ — frot of a source located at (n/9, 7r/4) with a resolution of 9.67 x 1CT 8 . 
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The FT of FM and complete response of the detector for a data of 120 days 
and 

f a = 25 Hz, h = h x = l \ 

a = 7 r/6, /3 = 7r/3, 7 = 27r/3, > (4.30) 

9 = n/9, = vr/4, ^ = tt/4. J 

are plotted in Figs. and (|4~T7T|l with a resolution of 1 /2T w 9.67x l()- s Hz. 
The corresponding power spectra are plotted respectively in Figs. (j4.11j) and ()4.12j) . 
Figs. ()4.13j) . ()4.14j) . (j4.15j) and ()4.16|) represent the power spectra of complete 
response of the detector at / + 2f rot , f — 2f rot , f + f ro t and / — f rot whereas 
Fig. f)4. 17j) shows the response of the detector at /. In this case also the most 
of the power of the signal lies in the frequency / + 2f rot and the least with 
/ — frot- The codes for computing the frequency modulated signal and com- 
plete response of the detector, written in FORTRAN, are given in appendices 
A and B. 

4.4 Spin down 

Pulsars loose their rotational energy by processes like electro-magnetic break- 
ing, emission of particles and emission of GW. Thus, the rotational frequency 
is not completely stable, but varies over a time scale which is of the order 
of the age of the pulsar. Typically, younger pulsars have the largest spin 
down rates. Current observations suggest that spin down is primarily due to 
electro- magnetic breaking (Manchester, 1992 and Kulkarni, 1992). Over the 
entire observing time T the frequency drift would be small but it may be 
taken into account for better sensitivity. To account this aspect we consider 
the evaluation of FT in a sequence of time windows by splitting the interval 
— T Q in M equal parts, each of interval At (T a = M A t) such that the 
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signal over a window may be treated as monochromatic. The strategy is to 
evaluate the FT over the window and finally to add the result. This process 
has been suggested by Brady and Creighton (2000) and Schutz (1998) in nu- 
merical computing and has been termed as stacking and tracking. For any 
such window let the time interval of data under consideration corresponds to 
t — t Q + n At and t — t a + (n + 1) At where t a is the instant representing the 
start of the data set and < n < M — 1. The window under consideration is 
the n t h window. Now we have 

et +(n+l)At 



/ h{i)e- i2nfI di 

Jt +nAt 
•At 



= / h(t + t + n A t)e- t27Tf ( t+t ° +nAt Ut ; 
Jo 



t = t + t + nAt 



(4.31) 
(4.32) 



Hence, the FT to account spin down is given via 



h(f) = / cos[$(t + t + n At)]e' i2n ^ t+to+nAt ^dt (4.33) 
J s Jo 



Taking the initial time of the data set 



t o = 



(4.34) 



and proceeding in the same manner as in previous section the FT is obtained 

as 



kf) 



i 



,i[27m(/ -/)At- 



2w rot 

k=oo 

2 J k{Z)i k cos k(a£- A) 



n - Q] e iv *[J (Z) + 
Jo 



k=i 



X 



J {N) + 2 J™Wi m cosm(e - C) 



m=l 



di (4.35) 
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After integration we get, 

Hf) 



V 



k=oo m=oo 



2w. 



rot 



E ^b[c s -w s 



k=—oo m=—oo 



where 



B 



X 

To 



(M^L + 2 nnAt(f -f)-n-Q 

J k (Z)J m {M) 

u 2 — (ak+m) 2 ' 

sin(z/r) cos(akr + mr — kX — m() 

-^{cos(^r) sin(afcr + mr — &A — raC) 

+ sin(/cA + mC)} , 

cos(z/t) cos(a/cr + mr — kX — m£) 
_l_afc±m s j n ^ r ^ s in(aArr 4- m r — kX — rn£) 

— cos(kX + ra£) , 
— anr , C = 5 — nr 
n 



r = £ = w rot T . 



-- w rot At, n = 0,1,2,3, ,M-1. 

The FT of the complete response would now be given via 



R(f) 



-i2/3 



(4.36) 



(4.37) 



h(f + 2f rot )/2 [h 0+ (F 1+ +iF 2+ )+h 0x (F 2x -iF lx )] + 

- s 

Hf - 2f rot )/2] [h D+ (F 1+ - iF 2+ ) - h 0x (F 2x + iF lx )] + 

- s 

h(f + f rot )/2\ [h 0+ (F 3+ + tF 4+ ) + h 0x (F Ax -iF 3 J] + 

- s 

h(f-frot)/2\ [h 0+ (F 3+ -zF 4+ )-h 0x (F, x +tF 3x )} + 



Hf) [h 0+ F 5+ -ih 0x F 5x ] 



(4.38) 



4.5 N-component signal 

The FT in Eqs. (J4.25J) and (|4.29|) are for a pulsar which emits GW signal at 
single frequency f Q . But there are known physical mechanisms which generate 
GW signals consisting of many components. An axially symmetric pulsar un- 
dergoing free precession, emits quadrupole GW at two frequencies, one equal 
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to the sum of the spin frequency and the precession frequency, and the other 
twice of it (Zimmermann, 1979, 1980). The quadrupole GW from a triaxial 
ellipsoid rotating about one of its principal axes consists of one component 
only (Thorne, 1987). In this case the signal has frequency about twice the 
spin frequency of the star. In general, if a star is non-axisymmetric and pre- 
cesses, the GW signal consists of more than two components. For the case 
of triaxial ellipsoid and small wobble angle there is a third component with 
frequency equal to twice the spin frequency of the star (Zimmermann, 1979, 
1980). Recently, new mechanisms e.g. r-mode instability of spinning neutron 
stars (Anderson, 1998; Lindblom et al., 1998; Owen, et al., 1998) and tem- 
perature asymmetry in the interior of the neutron star with miss-aligned spin 
axis (Bildsten, 1998) have been discussed in the literature. 

In view of the above discussion CGW signal may consists of frequencies 
which are multiple of some basic frequencies. An analysis of the GW data 
of N-component of signal has been made recently by Jaranowski and Krolak 
(2000). We in this section present Fourier analysis of N-component CGW 
signal. We model the N-component signal as 

N 

h{t) = 5>(0; (4-39) 
1=1 

h l+ (t) = /i 0j+ cos[$,(*)] (4.40) 
h lx (t) = h 0lx sin[^(t)] 1 = 1,2, ,N; (4.41) 

= W/ [t + Zi cos« rot - <j>) + Ni cos(£ rot - 5 t ) -Tlx- Qi] (4.42) 
where fi represent the component frequency of the signal. Let us write 

h l (t)=cos[$ l (t)] (4.43) 
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It is trivial matter to obtain hi(f), Ri(f) and to get 

M/) = £W) and (AAA) 
i 

R N (f) = ( 4 - 45 ) 

i 

We do not continue this analysis as the requisite formalism is analogous to what 
we have presented in Chapter 3 and earlier sections of the present Chapter. 



4.6 Discussion and conclusions 

The analysis and results obtained in earlier Chapter regarding FT of the re- 
sponse of a Laser Interferometer have been generalised in the present Chapter. 
In this context following points must be noted. 

1. For longer observation time, say, 120 days the resolution provided by 
FFT (equal to 1/T C ) is sufficient to represent the structure of side bands. 

2. Although in every case discussed, it turned out that the maximum power 
lies in the frequency / + 2f rot . However, this is not established conclu- 
sively whether this result is generic. In any case, for any given detector 
location, we can know in advance if similar type of behaviour is obtained. 

3. The computer time required in computing analytical FT depends on 
the orientation of the source and is independent of time interval of the 
observation data set. In contrast to this the computation cost increases 
for FFT with increase of the data set. 

4. Throughout our analysis in Chapters 3 and 4 we have employed following 
conditions. 

(i) The phase of the wave is zero at t — 0. 
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(ii) The observation time of the data set is from t — to t — T Q . 

5. As remarked in Chapter 3 the requirement 4-(i) may be acheived by 
translation of time-origin. Now let us see how the condition 4(h) can be 
relaxed. Let the data set is taken for 

t = ti to t = U + T Q (4.46) 

Eqs. fjPTJl . (IQ2l . and (PDMl) reveal that the results obtained 

there corresponds to data set 

t = n A t to t = (n + 1) A t (4.47) 

To obtain the sought generalisation we set 

ti = nAt and T Q = At (4.48) 

and may use the results of section ()4.3|) by noting that for the present 
case we have 

A = — nW or b A t =^> A = — Worbti (4.49) 

( = S - nw rot A t =^ ( = 8- w rot ti (4.50) 
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Chapter 5 

All sky search: Study of 
templates 

5.1 Introduction 

Gravitational wave Laser Interferometer antennas are essentially omni - di- 
rectional with their response better than 50% of the average over 75% of the 
whole sky (Grishchuk et al., 2000). Hence the data analysis systems will have 
to carry out all sky searches for the sources. We know that the amplitude of in- 
tense GW believed bathing Earth is very small, as compared to the sensitivity 
of GW detectors and is further masked by the dominant noise. In these cir- 
cumstances, CGW sources are of prime importance because for such sources 
we can achieve enhanced SNR by investigating longer observation data set. 
However, a long observation time introduces modulation effects, arising due 
to relative motion of detector and the source. As a consequence, there results 
distribution of power in a forest of sidebands to such an extent that there is 
reduction to the tune of 90% of the expected power due to AM. The problem 
of all sky search gains another dimension in view of the fact that there are rea- 
sons to believe the presence of intense GW sources whose locations and even 
frequencies are not known. Amongst such sources pulsars occupy an important 
position. Similar to all sky search one will also have to do all frequency search. 
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All sky all frequency search is the holy grail of gravitation pulsar astronomy. 
In this Chapter we confine ourselves to the problem of all sky search. 

Search of CGW without a priori knowledge appears to be computationally 
quite demanding even by the standard computers expected to be available 
in near future. For example, in the case of bandwidth 10 3 Hz, observation 
time 10 7 sec. and star's minimum decay time of 100 years one would require 
10 14 T flops computer (Frasca, 2000). Very fast computer and large memories 
with ample amount of disk space seems inevitable. However, choice of optimal 
data processing and clever programming is also integral part of solution to 
this problem. Amongst these the pre-correction of time series due to Doppler 
modulation before the data is processed may be a method, which will reduce 
computational requirements. In reference to this, Schutz (1991) has introduced 
the concept of patches in the sky defined as the region of space for which 
similar Doppler corrections would be required. He has also demonstrated that 
the number of patches required for 10 7 sec. observation data set and one 
KHz signal would be about 1.3 x 10 13 if one takes into account the rotational 
motion of Earth. However, the size of the patch would also depend on the 
data analysis technique being employed. 

Matched filtering is the most suitable technique for detection of signals 
from sources viz., pulsars whose wave form is known. The wave forms are 
used to construct a bank of templates, which represent expected signal wave 
form with all possible ranges of its parameters. The time of arrival, source 
location, frequency of the signal, ellipticity of the source and its spin down 
represent important parameters of GW emitted by a pulsar. For detection 
of GW we set a criterion to judge if the cross correlation of the templates 
with the corresponding data set exceeds the preassigned threshold. We would 
perform the data analysis employing the criterion of the FF. 



114 



5.2 Matched filter analysis: Templates 



The bank of templates will be matched to only a discrete set of signals from 
among the continuum of possible signals. Consequently, it is natural that all 
signals will not get detected with equal probability. However, it is possible 
to choose judiciously the set of templates so that all signals of a given ampli- 
tude are detected with a given minimum detection probability. The standard 
measure for deciding what class of wave form is good enough is the FF (Apos- 
tolatos, 1995). It quantitatively describes the closeness of the true signals to 
the template manifold in terms of the reduction of SNR due to cross corre- 
lation of a signal outside the manifold with all the templates lying inside the 
manifold. If the FF of a template family is unity the signal lies in the manifold. 
If the FF is less than unity the signal lies outside manifold. 

Even if the signal discrete templates lies within the template manifold it 
would be unlikely that any of the actual templates used would correspond to 
the signal. The parameters describing the search template (source location, 
ellipticity, etc.) can vary continuously through out a finite range of values. 
The set of templates characterised by the continuously varying parameters is 
ofcourse infinite. However, in practice the interferometer output must be cross 
correlated with a finite subset of the templates whose parameter values vary 
in discrete steps from one template to the next. This subset ("the discrete 
template family" ) has measure zero on the manifold of the full set of possible 
templates ("the continuous template family"), so the template which most 
closely match a signal will generally lie in between the signal and the nearest 
of the discrete template family. The mismatch between the signal and the 
nearest of the discrete templates will cause some reduction in SNR. This would 
mean that the members of the discrete template family must be chosen so as 
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to render acceptable loss of SNR. 

The study of templates has been made by many research workers in time 
domain [Schutz (1991), Krolak (1997), Brady et. al. (1998), Brady and 
Creighton (2000), Jaranowski et al. (1998) and Jaranowski and Krolak (2000)]. 
However, the analysis in frequency domain has the advantage of incorporat- 
ing interferometer's spectral noise density. In order to determine the number 
of templates required to perform matched filtering analysis let us rewrite the 
formula expressing FF [Eq. ()2.12|) ] as 



FF(e A) = W/)M/;fr.fr)> (5 . i; 

».« V(M/;0T,*r)IM/;0r.'M>W)IM/)> 



where, h(f) and h T (f; & T , 4> T ) represent respectively the FTs of the actual 
signal wave form and the templates. The template parameter 9t and 0t are 
chosen differing from the actual orientation of the source in discrete steps. 

We have seen in earlier Chapters that the AM of CGW data output results 
into redistribution of power at four additional frequencies / ± 2f rot , f ± f ro t in 
accordance with the FM. Hence it is sufficient for the analysis of FF to consider 
only the frequency modulated FT. The results obtained in Chapter 4 regarding 
FT of frequency modulated data output [Eq. (|4.25jl ] may be arranged using 
the symmetry property of the Bessel functions and one may write 



V 



W TO t 



J {Z)J {M) 



2u 2 



[{sin(ft + Q) - sin(ft +Q- uQ} + 



i {cos(ft + Q) - cos(ft + Q- v( )}] + 

J {z) V [(W - xv) -i{xu + yv)\ - 



m=l 
k=oo m=oo 

E E ^ 

k=l m=—oo 



C-iV 



(5.2) 
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X = sm{TZ + Q-mn/2) 

y = cos(K+ Q-mn/2) 

U = sin z/£ cos 77i(£ — 5) — ^ cos v^ Q sin m(£, 

V = cos u^ cos m(£ — 5) + ^7 sin z/£ sin m(£, 

Now it is straight forward to compute FF. To understand the procedure let us 
assume that there is a source at the location (9, 0) = (25°, 30°) emitting fre- 
quency f a = 0.5 Hz. We wish to analyze the data set for T Q = one sidereal day. 
We first fix the template parameter for the whole set given by 0t = = 30° 
and vary $t in discrete steps over its entire range i.e. 0° to 180°. The results 
obtained are plotted graphically in Fig. (|5.1|) . It is remarked that in order to 
compute the inner product defined via Eq. ()2.9j) one would require to integrate 
the expression over the band width of Doppler modulated signal. This may 
be determined either analytically by computing the maximum value of the 
Doppler shift in accordance with Eq. (|3.25J) or may be taken as represented by 
the frequency spectrum of the FT. In the present case we have taken the band 
width equal to 0.002 Hz. In a similar manner one may fix the parameter 
of the template set and obtain the variation of FF with template parameter 
4>t- Figures ()5.2|) and (|5.3J) represent respectively the behaviour for f Q = 25 
Hz, 6 = 9t = 1°, T D = one sidereal day for source orientation <p = 35° and 220°. 



— 5) — sin 77i(5 

— 5) — cos mS 
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The following points in reference to these plots may be noted. 

(i) The FF is unity for 9 T = 25°, 155° [Fig (JHU for T = 35°, 145° [Fig Q ] 
and for <p T = 220°, 320° [Fig IjOji ]. 

(ii) It is observed that decreases exponentially with template variables 
6^ and The specific relations are found numerically to be 

FF = e oM7(e T ~w (5 _ 4) 

FF = e°- 06( * T ~ ^ (5.5) 

(iii) The oscillatory behaviours do not represent any real situation as it arise 
because of the improper choice of numerical integration technique. How- 
ever, we are content with the technique we have employed as the region 
of such artificial facets fall in the region of the FF < 0.25. 

Finally, we conclude this section by noting the symmetry property of the tem- 
plate parameters. A closer look of the graphs and the remark (i) above reveal 
the following symmetry property. The FF is symmetrical under following 
t r ansf or mat ions . 

9 T — ► 7T - 9 T 0<9 T <n (5.6) 

(f>T ► K — <pT < (f)T < VT (5-7) 

(f)T ► 37T — <j)T 7T < (f)T < 271 (5-8) 

Let us note that these symmetry properties are based on our results for one 
sidereal day observation time. The generic nature of the symmetries may be 
established only after studying the variation of FF with T Q . Unfortunately, 
we could not make this analysis because of our limitations on computational 
facilities. 
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5.3 Number of Templates 

It is important to study the problem of number of templates for all sky search in 
the light of FF. The results of previous section reveal that the grid spacing A6 
in the parameter of templates may be expressed symbolically as a function 
of FF, f Q , and T Q i.e. 

A6 = F(FFJ ,T ) (5.9) 

Similarly, we have 

A<j) = g(FFJ ,T ) (5.10) 

In view of this, Equations ()5.4j) and (J5.5)) may be equivalently expressed as 

F(FF, 0.5, T) = (0.047) -1 ln(FF) (5.11) 

Q(FF, 25, T) = (0.06)- 1 ln(FF) (5.12) 

For any chosen value of FF one can determine A6 and A0. But there is 
no unique choice for it. Our interest would be in the assignment of AO and 
A0 such that the spacing are maximum resulting into the least number of 
templates. As we have mentioned earlier there is stringent requirement on 
reducing computer time. Accordingly, there is serious need of adopting some 
procedure/formalism to achieve this. For example, one may adopt the method 
of hierarchical search given by Mohanty and Dhurandhar (1996) and Mohanty 
(1998). This search is carried out in two steps. At the first level one would 
start with template bank with a coarse spacing in the parameter space but 
with a lower threshold. In the next level a more finely spaced set of templates 
and a higher threshold would be used but only around those templates of the 
previous level which crossed the previous threshold. 

However, an important issue related to the problem of number of templates 
is regarding the study of the behaviour of number of templates with FF for 
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different f a and T . We have made investigations of this aspect. We assume 
a source location (9, (ft) = (1°, 30°). We choose some value of FF, say 0.995. 
Now taking = 30° we determine the spacing A6 to yield the selected FF. 
In the case under investigation A8 is found to equal 4.5 x 10~ 5 . Thereafter, we 
introduce spacing A0 in the so obtained bank of templates and determine the 
resulting FF. The results obtained may be expressed in the form of a graph 
such as shown in Figs. (|5.4|) and (|5.5)1 . Interestingly the nature of these curves 
are similar. We have obtained a best fit to the graphs and obtain the relation. 

N T em P iates = exp[a - bx + ex 2 - dx 3 + ex 4 ] ; x = FF , (5.13) 

0.85 < x < 0.99 

where a, b, c, d and e are constants. The values of these constants are given 
in Table (lo~TJ) . 

Let us note from the graphs, for sake of comparison, that the number of 
templates required for FF equal to 0.97 are respectively 1.44 x 10 10 , 3.5 x 10 10 
and 5.5 x 10 10 for observation data set of 30, 120 and 365 days and f Q = 50 Hz. 
Similarly the number of templates required to analyse the observation data set 
of 120 days of GW frequencies 20, 50 and 100 Hz are respectively 1.22 x 10 10 , 
2.16 x 10 10 and 5 x 10 10 . It is observed that higher FF requires exponentially 
increasing large number of templates. 
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fo 

(Hz) 


T 

A o 

(days) 


a 

xl0~ 2 


b 

xlO- 4 


c 

xlO" 5 


d 

xl0~ 6 


e 

xl0~ 7 




30 


2138.05 


2071.43 


7225.73 


6239.43 


2036.14 


50 


180 


2317.05 


-71.3155 


1746.61 


2146.55 


944.931 




365 


2382.96 


216.917 


2464.42 


2464.42 


1031.54 


20 




2047.55 


-794.473 


3564.56 


4650.68 


1945.87 


50 


120 


2266.59 


4269.44 


15655.0 


17509.5 


6484.51 


100 




2360.23 


-206.906 


1158.27 


1491.01 


733.520 



Table 5.1: Coefficients of the best fit graphs obtained for the number of tem- 
plates. 
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Figure 5.4: Variation of number of templates with FF for fixed f a at different 
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Figure 5.5: Variation of number of templates with FF for fixed T Q and of 
frequencies. 
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5.4 Discussion 



In view of the complexity of the FT which contains exponential as well as 
Bessel functions; one has to be careful in computing FF. We have found 
useful to employ the Romberg integration using Pade approximation. We 
have used (i) QROMO of numerical recipes instead of QROMB as the former 
takes care of singularities, and (ii) RATINT routine for Pade approximation. 
The basic structure of the code is given in the appendix A and C (FM.F & 
TEM.F). 

To compute FF the code has been written in FORTRAN. The Bessel func- 
tions are computed using the routine BESSJO(X), BESSJl(X) and BESSJ(N,X) 
of numerical recipes (Press et al. 1986). The BESS JO (X) and BESSJl(X) com- 
pute the Bessel function of first kind of zeroth and first order respectively for 
any value of X, while BESSJ(N,X) computes the Bessel function of first kind 
of arbitrary integeral order N for any value of X. 

We have noticed marked symmetries in all sky search in both 9 and (f> space 
for one day observation time. It has been found that any FF corresponds 
to two values each in 9t and (fix- Accordingly computation burden will be 
reduced by a factor of four. However, it is not clear whether the symmetry 
property can be established analytically as well. The source location, because 
of the symmetry, is uncertain. Some other analysis is to be adopted for getting 
the exact location. 

We have computed the number of templates required for all sky search for 
matched filtering assuming the noise power spectral density S n (f) to be flat. 
However, for realistic situation, the effect of the noise have to be taken into 
consideration. 

The issues of optimum template parameterization and placement, and the 
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related computational burden have been discussed in literature by several au- 
thors notably by Sathyaprakash and Dhurandhar (1991), Dhurandhar and 
Sathyaprakash (1994), Owen (1996), Apostolatos (1995, 1996), Mohanty and 
Dhurandhar (1996), Mohanty (1998), Owen and Sathayaprakash (1999). The 
question of possible efficient interpolated representation of the correlators is a 
problem of current interest and remains still unsolved. 
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Appendix A 



Source code of the frequency 
modulated signal: FM.F 

C FT of the frequency modulated signal 

FUNCTION hfm(fnu,dtheta,dphi) 
IMPLICIT none 
INTEGER k,m,ms,ks 
REAL t,fO,fnu,dtheta,dphi,dalpha,dbetaO 
C [t = T , fO = f a , fnu = v, dtheta = 9°, dphi = <p°, dalpha = a°. 

* dbetaO = (3°\ 

PARAMETER (t = 86164.0, fO = 25.0, dalpha = 75.0, 
dbetaO = 60.0, ks = 720, ms = 8) 
C ['ks' & 'ms' are the values of 'k' & 'm' respectively] 

REAL h,z,bessjO,bessjl,bessj,bsh(-ms:ms),bsz(0:ks) 
REAL*8 pi,fr,pfr,a,u,rc,rad,theta,phi,epsilon,alpha,betaO, 

* sth,sal,sbeO,sph,sep,cth,cal,cbeO,cph,cep,p,q,r,delta,crq, 

* srq,bs0,fnn,cfnn,sfnn,rqf,fnu2,rplusq 

REAL*8 dm(l:ms),dmu(l:ms),rqm(l:ms),cmud(l:ms), smud(l:n 

* smd(l:ms),cmd(l:ms),crqm(l:ms),srqm(l:ms),wmn(l:ms),vmn(l 

* akd(l:ks,-ms:ms),rqk(l:ks,-ms:ms),skm(l:ks,-ms:ms), 

* pmd(l:ks,-ms:ms),akm(l:ks,-ms:ms),ckmu(l:ks,-ms:ms), 

* skmu(l:ks,-ms:ms),ckmr(l:ks,-ms:ms),ckm(l:ks,-ms:ms), 

* skmr(l:ks,-ms:ms) 
COMPLEX suml,sum2,hfm 

C [hfm — > Spectrum of the frequency modulated signal.] 
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EXTERNAL bessjO,bessjl,bessj 



pi = 4.d0*datan(l.d0) 

fr = I.d0/86164.d0 

pfr = 2.d0*pi*fr 

a = I.d0/365.2572d0 

u = 2.d0*pi*fr*t 

rc = 2.d0*pi*f0*0.0637/3.d0 

rad = pi/180.d0 

theta = rad*dtheta 

phi = rad*dphi 

epsilon = rad*23.45d0 

alpha = rad* 75. dO 

betaO = rad*60.d0 

sth = dsin(theta) 

sal = dsin(alpha) 

sbeO = dsin(betaO) 

sph = dsin(phi) 

sep = dsin(epsilon) 

cth = dcos(theta) 

cal = dcos(alpha) 

cbeO = dcos(betaO) 

cph = dcos(phi) 

cep = dcos(epsilon) 

p = rc*sal*(cbeO*(sth*cep*sph + cth*sep) - sbeO*sth*cph) 

q = rc*sal*(sbeO*(sth*cep*sph + cth*sep) + cbeO*sth*cph) 

z = 2.d0*pi*f0*1496.d0/3.d0*sth 

r = z*cph 

rplusq = r + q 

delta = datan(p/q) 

h = dsqrt(p*p + q*q) 

crq = dcos (rplusq) 

srq = dsin(rplusq) 

bsz(O) = bessjO(z) 
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bsh(O) = bessjO(h) 
bsz(l) = bessjl(z) 
bsh(l) = bessjl(h) 
bsh(-l) = - bsh(l) 
bsO = bsz(0)*bsh(0) 

DO 10 m = 2,ms 
bsh(m) = bessj(m,h,bsh(0),bsh(l)) 
bsh(-m) = (-l)**m*bsh(m) 
CONTINUE 

DO 20 k = 2,ks 
bsz(k) = bessj(k,z,bsz(0),bsz(l)) 
CONTINUE 

DO 30 m = l,ms 
dm(m) = m*delta 
dmu(m) = m*(u-delta) 
rqm(m) = rplusq - m*pi/2.d0 
cmud(m) = dcos(dmu(m)) 
smud(m) = m*dsin(dmu(m)) 
smd(m) = m*dsin(dm(m)) 
cmd(m) = dcos(dm(m)) 
crqm(m) = dcos(rqm(m)) 
srqm(m) = dsin(rqm(m)) 
CONTINUE 

DO 40 m = -ms,ms 
DO 50 k = l,ks 
akd(k,m) = a*k*u+m*u-k*phi - m*delta 
rqk(k,m) = (k + m)*pi/2.d0 - rplusq 
pmd(k,m) = k*phi+m*delta 
akm(k,m) = k*a + m 
ckmu(k,m) = dcos(akd(k,m)) 
skm(k,m) = akm(k,m)*dsin(pmd(k,m)) 
skmu(k,m) = akm(k,m)*dsin(akd(k,m)) 
ckmr(k,m) = dcos(rqk(k,m)) 
ckm(k,m) = dcos(pmd(k,m)) 
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skmr(k,m) = dsin(rqk(k,m)) 
50 CONTINUE 
40 CONTINUE 

fnn = fnu*u 

fnu2 = fnu*fnu 

cfnn = dcos(fnn) 

sfnn = dsin(fnn) 

rqf = rplusq - fnn 

suml = cmplx(0. 0,0.0) 

sum2 = cmplx(0.0,0.0) 
DO 60 m = l,ms 

wmn(m) = sfnn*cmud(m) - (cfnn*smud(m) + smd(m))/fnu 

vmn(m) = cfnn*cmud(m) + sfnn*smud(m)/fnu - cmd(m) 

suml = suml + bsz(0)*bsh(m)/(fnu2 - m*m)* 

* cmplx((wmn(m)*crqm(m) - vmn(m)*srqm(m)), 
*- (wmn(m)*srqm(m) + vmn(m)*crqm(m))) 

60 CONTINUE 

DO 70 m = - ms,ms 
DO 80 k = l,ks 
sum2 = sum2 + bsz(k)*bsh(m)/(fnu2 - (a*k+m)**2)* 

* cmplx((sfnn*ckmu(k,m) - (skm(k,m) + skmu(k,m)*cfnn)/fnu), 
*- (cfnn*ckmu(k,m) - ckm(k,m) + skmu(k,m)* 

* sfnn / fnu) ) * cmplx (ckmr (k,m) , skmr (k, m) ) 
80 CONTINUE 

70 CONTINUE 

hfm = fnu/pfr*(bs0/(2.0*fnu2)*cmplx((srq - dsin(rqf)), 

* (crq - dcos(rqf))) + suml + sum2) 

RETURN 
END 
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Appendix B 



Source code of the noise free 
complete response of the 
detector: CR.F 

C FT of the nosie free complete response of the detector 

IMPLICIT none 

INTEGER k,m,n,ms,ks,j,l 

PARAMETER (ks = 27225,ms = 15,j = 174) 
C [ j — ► Number of frequencies bins in the spectrum 
C of the complete response.] 

REAL t,fO,f,dtheta,dphi,dalpha,dbetaO,dgamma,dpsi,hO_p, 

* hO_c,res_n,z,h,bessjO,bessjl,bessj 

PARAMETER (t = 86164.0, fO = 50.0, f = 49.9995, dtheta = 10.0, 
*dphi = 0.0, dalpha = 46.45, dbetaO = 0.0, dgamma = 171.8, 
*dpsi = 45.0, h0_p = 1.0, h0_c = 1.0, res_n = 2.0*t) 
C [dgamma = 7 , dpsi = h0_p & h0_c — > Amplitude of the plus 

C and cross polarisation of the signal respectively and 
C res_n — > resolution of the spectrum.] 

REAL*8 pi, rad,fr,pfr,a,u,rc,theta,phi,epsilon, alpha, 

* psi,betaO,sth,sal,sbeO,sph,sep,cth,cal,cbeO,cph,cep,gamma,ssi,csi 

* ,s2ep,c2ep,s2gm,c2gm,al,am,an,ap,aq,ar,ag,ah,ai,aj,fl_p,f2_p,f3_p 

* ,f4_p,f5_p,sph_c,cph_c,al_c,am_c,an_c,ap_c,aq_c,ar_c,ag_c,ah_c, 

* ai_c,aj_c,fl_c,f2_c,f3_c,f4_c,f5_c,p,q,r,rplusq,delta,crq,srq,fn(0:j) 
C [fn — > frequency bins of the complete response.] 
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REAL*8 bsh(-ms:ms),bsz(0:ks),bsO,dm(l:ms),dmu(l:ms), 

* rqm(l:ms),cmud(l:ms),smud(l:ms),smd(l:ms),cmd(l:ms),crqm(l:ms), 

* srqm(l:ms),akd(l:ks,-ms:ms),rqk(l:ks,-ms:ms),pmd(l:ks,-ms:ms), 

* akm(l:ks,-ms:ms),ckmu(l:ks,-ms:ms),skm(l:ks,-ms:ms), 

* skmu(l:ks,-ms:ms),ckmr(l:ks,-ms:ms), ckm(l:ks,-ms:ms), 

*skmr(l:ks,-ms:ms),fnu(-2:2,0:j),fnn(-2:2,0:j), 

*cfnn(-2:2,0:j),sfnn(-2:2,0:j),rqf(-2:2,0:j),wmn(l:ms,-2:2), 

* vmn(l:ms,-2:2) 

COMPLEX fp2fr,fm2fr,fpfr,fmfr,suml(-2:2,0:j),sum2(-2:2,0:j), 
*hfm(-2:2,0:j),crp2f(0:j),crplf(0:j),cr0f(0:j),crmlf(0:j), 
*crm2f(0:j),cr(0:j) 
C [cr — > complete response of the detector and crp2f, crplf, crOf, 
C crm2f & crmlf — > complete response at / + 2f r , f + f, n /, 
C / - 2/ r & / - f r respectively.] 
EXTERNAL bessjO,bessjl,bessj 

pi = 4.d0*datan(l.d0) 

rad = pi/180.d0 

fr = I.d0/86164.d0 

pfr = 2.d0*pi*fr 

a = I.d0/365.2572d0 

u = 2.d0*pi*fr*t 

rc = 2.d0*pi*f0*0.0637d0/3.d0 

theta = rad*dtheta 

phi = rad*dphi 

epsilon = rad*23.45d0 

alpha = rad*dalpha 

psi = rad*dpsi 

betaO = rad*dbetaO 

sth = dsin(theta) 

sal = dsin(alpha) 

sbeO = dsin(betaO) 

sph = dsin(phi) 

sep = dsin(epsilon) 
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cth = dcos(theta) 

cal = dcos(alpha) 

cbeO = dcos(betaO) 

cph = dcos(phi) 

cep = dcos(epsilon) 

gamma = rad*dgamma 

ssi = dsin(psi) 

csi = dcos(psi) 

s2ep = dsin(2.d0*epsilon) 

c2ep = dcos(2.d0*epsilon) 

s2gm = dsin(2.d0*gamma) 

c2gm = dcos(2.d0*gamma) 

Beam pattern for plus polarisation 

al = csi*cph - cth*sph*ssi 
am = csi*sph + cth*cph*ssi 
an = - ssi*cph - cth*sph*csi 
ap = - ssi*sph + cth*cph*csi 
aq = sth*sph 
ar = sth*cph 

ag = ((al*aq + am*ar)*sep - (al*an - am*ap)*cep)/2.d0 
ah = ((an*an - ap*ap)*cep*cep - (al*al - am*am) + 

* (aq*aq - ar*ar)*sep*sep - (an*aq + ap*ar)*s2ep)/2.d0 
ai = ((aq*aq - ar*ar)*(s2ep) - (an*an - ap*ap)*s2ep - 

*2.d0*(an*aq + ap*ar)*c2ep)/2.d0 
aj = ((al*aq + am*ar)*cep + (al*an - am*ap)*sep)/2.d0 
fl_p = h0_p*(ah*s2gm*(cal*cal + 1) - 4.d0*ag*cal*c2gm)/2.d0 
f2_p = h0_p*(ah*cal*c2gm + ag*s2gm*((cal)**2 + l.dO)) 
f3_p = h0_p*(ai*sal*c2gm + aj*dsin(2.d0*alpha)*s2gm) 
f4_p = h0_p*(4.d0*aj*sal*c2gm - ai*dsin(2.d0*alpha)*s2gm)/2.d0 
f5_p = h0_p*sal*sal*s2gm*(ah + al*al - am*am)*3.d0/2.d0 

Beam pattern for cross polarisation 
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sph_c = dsin(phi - pi/4.d0) 
cph_c = dcos(phi - pi/4.dO) 
al_c = csi*cph_c - cth*sph_c*ssi 
am_c = csi*sph_c + cth*cph_c*ssi 
an_c = - ssi*cph_c - cth*sph_c*csi 
ap_c = -ssi*sph_c + cth*cph_c*csi 
aq_c = sth*sph_c 
ar_c = sth*cph_c 

ag_c = ((al_c*aq_c + am_c*ar_c)*sep - (al_c*an_c - 

* am_c*ap_c)*cep)/2.d0 

ah_c = ((an_c*an_c - ap_c*ap_c)*cep*cep - 

* (al_c*al_c - am_c*am_c) +(aq_c*aq_c - ar_c*ar_c)*sep*sep - 

* (an_c*aq_c + ap_c*ar_c)*s2ep)/2.d0 

ai_c = ((aq_c*aq_c - ar_c*ar_c)*(s2ep) - (an_c*an_c - ap_c*ap_c)* 
*s2ep - 2.d0*(an_c*aq_c + ap_c*ar_c)*c2ep)/2.d0 
aj_c = ((al_c*aq_c + am_c*ar_c)*cep + (al_c*an_c - 

* am_c*ap_c)*sep)/2.d0 

fljc = h0_c*(ah_c*s2gm*(cal*cal + 1) - 4.d0*ag_c*cal*c2gm)/2.d0 
f2_c = h0_c*(ah_c*cal*c2gm+ag_c*s2gm*((cal)**2 + l.dO)) 
f3_c = h0_c*(ai_c*sal*c2gm + aj_c*dsin(2.d0*alpha)*s2gm) 
f4_c = h0_c*(4.d0*aj_c*sal*c2gm - 
*ai_c*dsin(2.d0*alpha)*s2gm)/2.d0 
f5_c = h0_c*sal*sal*s2gm*(ah_c + al_c*al_c - am_c*am_c)*3.d0/2.d0 

C Both polarisations are suitably arrange to compute 
C complete response 

fp2fr = cmplx((fl_p + f2_c), (f2_p - fl_c))/2.d0 
fm2fr = cmplx((fl_p - f2_c), (- f2_p - fl_c))/2.d0 
fpfr = cmplx((f3_p + f4_c), (f4_p - f3_c))/2.d0 
fmfr = cmplx((f3_p - f4_c), (- f4_p - f3_c))/2.d0 

C Computation for frequncy modulation 
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p = rc*sal*(cbeO*(sth*cep*sph + cth*sep) - sbeO*ar) 

q = rc*sal*(sbeO*(sth*cep*sph + cth*sep) + cbeO*ar) 

z = 2.d0*pi*f0*1496.d0/3.d0*sth 

r = z*cph 

rplusq = r + q 

delta = datan(p/q) 

h = dsqrt(p*p + q*q) 

crq = dcos(rplusq) 

srq = dsin (rplusq) 

bsz(O) = bessjO(z) 

bsh(O) = bessjO(h) 

bsz(l) = bessjl(z) 

bsh(l) = bessjl(h) 

bsh(-l) = - bsh(l) 

bsO = bsz(0)*bsh(0) 
DO 10 m = 2,ms 

bsh(m) = bessj(m,h,bsh(0),bsh(l)) 

bsh(-m) = (-l)**m*bsh(m) 
10 CONTINUE 

DO 20 k = 2,ks 

bsz(k) = bessj(k,z,bsz(0),bsz(l)) 
20 CONTINUE 

DO 30 m = l,ms 

dm(m) = m*delta 

dmu(m) = m*(u-delta) 

rqm(m) = rplusq - m*pi/2.d0 

cmud(m) = dcos(dmu(m)) 

smud(m) = m*dsin(dmu(m)) 

smd(m) = m*dsin(dm(m)) 

cmd(m) = dcos(dm(m)) 

crqm(m) = dcos(rqm(m)) 

srqm(m) = dsin(rqm(m)) 
30 CONTINUE 
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DO 40 m = - ms,ms 
DO 50 k = l,ks 
akd(k,m) = a*k*u+m*u - k*phi - m* delta 
rqk(k,m) = (k + m)*pi/2.d0 - rplusq 
pmd(k,m) = k*phi + m*delta 
akm(k,m) = k*a + m 
ckmu(k,m) = dcos(akd(k,m)) 
skm(k,m) = akm(k,m)*dsin(pmd(k,m)) 
skmu(k,m) = akm(k,m)*dsin(akd(k,m)) 
ckmr(k,m) = dcos(rqk(k,m)) 
ckm(k,m) = dcos(pmd(k,m)) 
skmr(k,m) = dsin(rqk(k,m)) 
50 CONTINUE 
40 CONTINUE 

DO 60 n = 0,j 
fn(n) = f + n/res_n 
60 CONTINUE 

DO 70 n = 0,j 
DO 80 1 = - 2,2 
fnu(l,n) = (fO - fn(n))/fr + 1 
fnn(l,n) = fnu(l,n)*u 
cfnn(l,n) = dcos(fnn(l,n)) 
sfnn(l,n) = dsin(fnn(l,n)) 
rqf(l,n) = rplusq - fnn(l,n) 
suml(l,n) = cmplx(0. 0,0.0) 
sum2(l,n) = cmplx(0.0,0.0) 

DO 90 m=l,ms 
wmn(m,l) = sfnn(l,n)*cmud(m) - (cfnn(l,n)*smud(m) + 

* smd(m))/fnu(l,n) 

vmn(m,l) = cfnn(l,n)*cmud(m) + sfnn(l,n)* 

* smud(m)/fnu(l,n) - cmd(m) 

suml(l,n) = suml(l,n) + bsz(0)*bsh(m)/(fnu(l,n)**2 - m*m)* 

* cmplx((wmn(m,l)*crqm(m) - vmn(m,l)*srqm(m)), - 

* (wmn(m,l)*srqm(m) + vmn(m,l)*crqm(m))) 
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90 CONTINUE 

DO 100 m = - ms,ms 
DO 110 k = l,ks 
sum2(l,n) = sum2(l,n) + bsz(k)*bsh(m)/(fhu(l,n)**2 - 

* (a*k + m)**2)* 

* cmplx((sfnn(l,n)*ckmu(k,m) - (skm(k,m) + skmu(k,m)*cfnn(l,n))/ 
*fnu(l,n)), - (cfnn(l,n)*ckmu(k,m) - ckm(k,m)+skmu(k,m)* 

* sfnn(l,n)/fnu(l,n)))*cmplx(ckmr(k,m),skmr(k,m)) 
110 CONTINUE 

100 CONTINUE 

hfm(l,n) = fnu(l,n)/pfr*(bs0/(2.d0*fnu(l,n)**2)* 

*cmplx((srq - dsin(rqf(l,n))),(crq - dcos(rqf(l,n)))) + 

*suml(l,n) + sum2(l,n)) 
80 CONTINUE 
70 CONTINUE 

C Computation for complete response 

DO 120 n=0j 
crp2f(n) = fp2fr*hfm(2,n) 
crplf(n) = fpfr*hfm(l,n) 
cr0f(n) = (cmplx(f5_p,f5_c))*hfm(0,n) 
crmlf(n) = fm2fr*hfm(-l,n) 
crm2f(n) = fmfr*hfm(-2,n) 

cr(n) = crp2f(n) + crplf(n) + cr0f(n) + crmlf(n) + crm2f(n) 
write(10,*) fn(n),cabs(cr(n)),cabs(crp2f(n))),cabs(crplf(n)) 
write(20,*) fn(n),cabs(crm2f(n)),cabs(crmlf(n))),cabs(cr0f(n)) 
120 CONTINUE 
STOP 
END 
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Appendix C 



Source code for the estimation 
of templates for all sky search: 

TEM.F 



C Estimation of templates for all sky search 

IMPLICIT none 
INTEGER n,nt 
PARAMETER (nt = 721) 
REAL*8 fr,pi 

REAL x,y,fO,st,sig,dtheta,dphi,tdtheta,tdphi(l:nt),sigsig,ff(l:nt), 

* sigtem(l:nt),temtem(l:nt),tem_theta,tem_phi(l:nt),skytem(l:nt) 

C [ff = Fitting Factor, skytem = Number of templates for all sky search] 

C FREQUENCY AND SOURCE LOCATION 

PARAMETER (fO = 1.0, dtheta = 1.0, tdtheta = 0.999955, 

* dphi = 35.0) 
EXTERNAL st,qromo,sig 

pi = 4.d0*datan(l.d0) 
fr = I.d0/86164.d0 
x = (fO - 0.999998d0)/fr 
y = (fO - 1.000002d0)/fr 

CALL qromo(sig,x,y, dtheta, dphi, dtheta,dphi,sigsig) 
tem_theta = 180.0/(dtheta - tdtheta) 
DO 10 n = 0,nt 
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tdphi(n) = 0.5*float(n) 

tem_phi(n) = abs(360.d0/(dphi + 0.1e-4 - tdphi(n))) 
CALL qromo(sig,x,y,tdtheta,tdphi(n),tdtheta,tdphi(n),temtem(n)) 
CALL qromo(st,x,y,dtheta,dphi,tdtheta,tdphi(n),sigtem(n)) 
ff(n) = abs(sigtem(n)/(sqrt(sigsig*temtem(n)))) 
skytem(n) = tem_theta*tem_phi(n) 
write(lO) sngl(ff(n)),skytem(n) 
10 CONTINUE 
STOP 
END 

SUBROUTINE qromo(func,a,b,dtheta,dphi,tdtheta,tdphi,ss) 
INTEGER JMAX,JMAXP,K,KM 
REAL a,b,func,ss,EPS,dtheta,dphi,tdtheta,tdphi 
EXTERNAL func 

PARAMETER (EPS = l.e-6, JMAX = 14, JMAXP = JMAX+1, 
*K = 5, KM = K-l) 
C USES ratint,midpnt 
INTEGER j 

REAL dss,h( JMAXP) ,s( JMAXP) 
h(l) = 1. 

DO 10 j = 1,JMAX 
CALL midpnt(func,a,b,dtheta,dphi,tdtheta,tdphi,s(j),j) 

IF (j.ge.K) THEN 
CALL ratint(h(j - KM),s(j - KM),K,0.,ss,dss) 

IF (abs(dss)Je.EPS*abs(ss)) RETURN 

ENDIF 
s(j+l) = s(j) 
h(j + l) = h(j)/9. 
10 CONTINUE 

PAUSE 'too many steps in qromo' 

END 
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SUBROUTINE midpnt (func,a,b,dtheta,dphi,tdtheta,tdphi,s 
INTEGER n 

REAL a,b,s,func,dtheta,dphi,tdtheta,tdphi 

EXTERNAL func 

INTEGER it J 

REAL ddel,del,sum,tnm,x 

IF (n.eq.l) THEN 
s = (b - a)*func(0.5*(a + b),dtheta,dphi,tdtheta,tdphi) 

ELSE 
it = 3**(n - 2) 
tnm = it 

del = (b - a)/(3.*tnm) 
ddel = del + del 
x = a + 0.5*del 
sum = 0. 

DO 10 j=l,it 
sum = sum + func(x,dtheta,dphi,tdtheta,tdphi) 
x = x + ddel 

sum = sum + func(x,dtheta,dphi,tdtheta,tdplii) 
x = x + del 
CONTINUE 

s = (s + (b - a)*sum/tnm)/3. 
ENDIF 
RETURN 
END 



SUBROUTINE ratint(xa,ya,n,x,y,dy) 
INTEGER n,NMAX 
REAL dy,x,y,xa(n),ya(n),TINY 
PARAMETER (NMAX = 10,TINY = l.e-25) 
INTEGER i,m,ns 

REAL dd,h,hh,t,w,c(NMAX),d(NMAX) 
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ns = 1 

hh = abs(x - xa(l)) 

DO 10 i = l,n 
h = abs(x - xa(i)) 

IF (h.eq.O.) THEN 

y = ya(i) 

dy = 0.0 
RETURN 

ELSE IF (h.lt.hh) THEN 
ns = i 
hh = h 

ENDIF 
c(i) = ya(i) 
d(i) = ya(i) + TINY 
10 CONTINUE 
y = ya(ns) 
ns = ns - 1 

DO 20 m = l,n - 1 
DO 30 i = l,n - m 
w = c(i + 1) - d(i) 
h = xa(i + m) - x 
t = (xa(i) - x)*d(i)/h 
dd = t - c(i + 1) 

IF(dd.eq.0.)PAUSE 'failure in ratint' 
dd = w/dd 
d(i) = c(i + l)*dd 
c(i) = t*dd 
30 CONTINUE 

IF (2*ns.lt.n - m) THEN 
dy = c(ns + 1) 

ELSE 
dy = d(ns) 
ns = ns - 1 

ENDIF 



144 



y = y + dy 

CONTINUE 
RETURN 
END 

FUNCTION sig(fnu,dtheta,dphi) 
REAL sig,fnu,dtheta,dphi 
COMPLEX hfm 
EXTERNAL hfm 

sig = real(hfm(fnu,dtheta,dphi)*conjg(hfm(fnu,dtheta,dphi))) 
RETURN 
END 

FUNCTION st(fnu,dtheta,dphi,tdtheta,tdphi) 
REAL st,fnu,dtheta,dphi,tdtheta,tdphi 
COMPLEX hfm 
EXTERNAL hfm 

st = real(hfm(fnu,dtheta,dphi)*conjg(hfm(fnu,tdtheta,tdphi))) 
RETURN 
END 
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